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SUMMARY 


Significant  improvements  in  the  capacity  and  service  life  of  reinforced 
earth  structures  require  an  improved  understanding  of  the  fundamental  behavior 
of  these  systems.  Both  experimental  and  analytical  investigations  were  carried 
out  to  develop  models  for  the  interaction  of  geotextile-type  reinforcement  and 
granular  soils.  Reinforcement  configurations  and  systems  investigated  were 
thought  to  be  applicable  to  temporary  runways,  also  known  as  alternate  launch 
and  recovery  surfaces  (ALRS).  Model  ALRS  systems  using  geotextiles  and  geogrids 
as  reinforcement  were  tested  in  the  laboratory  in  a  variety  of  configurations. 

These  were  loaded  to  failure,  quasi-statically,  by  both  plane  strain  and  axisym- 
metric  rigid  plates.  Load-deformation  characteristics,  as  well  as  the  shape  of 
the  deflected  basin,  are  reported. 

Two  new  analytical  models  for  fabri c-rei nforced  soil  systems  were 
developed.  The  first  is  based  on  probabilistic  concepts  for  vertical  stress 
diffusion  in  a  particulate  medium.  Numerical  solutions  were  obtained  for  some 
practical  subgrade-rei nforci ny  situations.  The  efficiency  of  the  membrane  was 
found  to  increase  as  the  underlying  subgrade  becomes  more  compressible,  as  long 
as  the  geotextile  is  sufficiently  strong  and  possesses  sufficient  frictional 
resistance.  Another  model  was  developed  which  considers  the  diffusion  of  strain 
through  a  particulate  medium.  This  model  predicts  surface  deflection  profiles 
as  well  as  strains  with  depth,  and  the  approach  has  considerable  potential  in 
foundation  engineering  as  well  as  in  soil -fabric  reinforcement  systems. 

Significant  increases  in  bearing  capacity  and  modulus  of  subgrade  reaction 
as  a  function  of  both  the  location  (depth)  of  the  reinforcement  and  the  number 
of  reinforcement  layers  were  observed.  However,  there  was  a  decrease  in 
improvement  as  the  depth  to  the  first  layer  increased.  The  "critical  depth,"  if 
it  exists,  is  probably  about  one-third  the  width  (or  diameter)  of  the  loaded  area. 
Edye  fixity  conditions  were  found  to  be  relatively  unimportant,  and  the  benefit 
of  multiple-reinforcement  layers  was  found  to  be  greater  if  the  depth  and  spac¬ 
ing  were  small  compared  to  the  diameter  of  the  loaded  area.  Surprisingly,  lit¬ 
tle  difference  in  response  between  geogrids  and  yeotextiles  was  observed,  prob¬ 
ably  because  the  sand  used  in  the  experiments  did  not  provide  the  interlock  com¬ 
ponent  apparently  necessary  for  the  optimum  functioning  of  the  geogrids.  Scaling 
of  maximum  load,  bearing  capacity,  etc.,  based  on  the  diameters  of  the  loaded 
areas  was  unsuccessful.  Surface  deflection  profiles  due  to  the  loaded  plate  at  one- 
half  maximum  load  could  be  reasonably  predicted  from  the  strain  diffusion 
hypothesis  developed  earlier. 

Implications  of  the  research  findings  for  practical  ALRS  systems,  as  well  as 
an  outline  of  plans  for  Phase  II  research,are  also  presented. 
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SECTION  I 


INTRODUCTION  AND  TECHNICAL  BACKGROUND 

Increased  understanding  of  the  fundamental  behavior  of  certain  mechanical 
reinforcing  and  stabilization  techniques  is  necessary  before  significant 
improvements  can  be  made  in  the  load-carrying  capacity  and  service  life  of  rein¬ 
forced  earth  structures.  Substantial  advancements  would  be  expected  in  several 
aspects  of  airfield  pavement  systems,  particularly  in  expedient  or  alternate 
launch  and  recovery  surfaces  and  for  rapid  repair  of  bomb-damaged  runways.  The 
knowledge  gained  is  also  appropriate  to  other  potential  application  areas  such 
as  retaining  walls,  embankments,  and  protective  construction. 

I.  EARTH  REINFORCING 

Soils  often  require  reinforcing  to  accommodate  anticipated  loadings.  ’'Cor¬ 
duroy"  roads,  constructed  on  logs  or  timbers  in  colonial  North  America  and  Scan¬ 
dinavia,  are  a  form  of  reinforced  earth,  as  are  bamboo  fascines  used  under  low 
embankments  in  southeast  Asia.  Embankments  have  often  been  constructed  directly 
on  the  brush  and  small  trees  which  are  common  on  marshy  lands.  In  recent  years, 
a  reinforcing  system  of  two  rows  of  short  sheet  piles  or  steel  channel  sections 
connected  by  steel  anchor  rods  has  been  developed  at  the  Swedish  Geotechnical 
Institute  to  increase  the  stability  of  embankments  constructed  on  soft  founda¬ 
tions  (Reference  1).  Reinforcing  has  also  been  carried  out  using  woven  and 
nonwoven  fabrics,  plastic  and  steel  nets,  used  automobile  tire  casings,  ordinary 
landing  mats,  "Columbus"  fascine  mats  (Reference  2),  or  reinforced  plastic  or 
rubber  membranes.  Reference  3  summarizes  many  of  these  recent  developments  for 
reinforcing  both  embankments  and  retaining  walls. 

The  primary  functional  requirements  of  the  reinforcement  elements  are  that: 
(1)  they  must  have  a  sufficiently  high  deformation  modulus  in  tension;  and  (2) 
they  must  be  able  to  develop  sufficient  frictional  resistance  with  the  subsoil 
and/or  embankment  materials.  Many  plastics  and  nonwoven  fabrics  have  creep  pro¬ 
perties  such  that  their  effectiveness  as  reinforcement  may  decrease  with  time. 

Of  course,  in  reinforced  embankments  the  strength  of  the  subsoil  may  increase 
faster  than  the  corresponding  creep  in  the  rei nforcement  to  neutralize  the 
effect  of  creep.  Or,  in  the  case  of  reinforced  sands  or  retaining  wall  back¬ 
fills  in  which  consol i dati on  does  not  take  place,  creep  deformations  under  high 
loads  could  be  significant.  However,  for  the  relatively  short-term  loadings 
considered  in  this  research,  the  creep  of  plastic  reinforcement  materials  was 
thought  to  present  no  particular  difficulty. 

The  use  of  woven  and  nonwoven  fabric  materials  (ASTM:  "geotextiles")  is  a 
relatively  recent  development  in  the  U.S.  With  a  few  notable  exceptions,  most 
of  the  research  and  development  work  on  these  materials  was  done  in  Europe,  and 
applications  were  primarily  directed  toward  stabilizing  temporary  roads  on  soft 
foundations.  During  the  past  10  years,  many  European  nonwoven  fabrics  have 
become  available  in  the  U.S.,  and  their  use  for  certain  specific  situations  is 
i ncreasi ng. 

With  only  a  few  exceptions,  woven  geotextile  technology  began  in  the  U.S. 
and  then  moved  to  Europe.  Initial  applications  in  the  U.S.  were  with  woven 
monofilament  fabrics  used  as  "filter  fabrics,"  tnat  is,  as  an  alternative  to 


granular  filters  under  riprap  and  in  other  erosion  control  features.  Recent 
developments  in  woven  technology  have  included  the  slit  film  fabrics,  which  are 
much  stronger  and  have  a  higher  modulus  than  typical  nonwovens,  but  cost  about 
the  same  per  unit  area. 

Considerable  research  using  woven  fabrics  as  reinforcing  has  been  conducted 
at  the  Swedish  Geotechnical  Institute  and  at  Purdue  University.  This  work 
involved  laboratory  investigations  (References  4  and  5)  and  three  field  instal¬ 
lations  using  a  woven  polyester  fabric  (References  6  and  7).  This  same  fabric 
was  used  in  the  present  investigations. 

Other  U.S.  research  on  woven  fabrics  has  been  carried  out  primarily  by  T. 

A.  Haliburton  and  his  associates.  They  were  interested  in  the  problem  of 
embankments  constructed  on  very  soft  foundations  (References  8,  9,  and  10).  They 
have  also  conducted  experimental  research  on  the  mechanical  properties  of 
geotextiles  (Reference  11)  and  on  soil-fabric  interaction  (Reference  12).  Of 
considerable  interest  to  the  present  research  is  Reference  13,  in  which  the 
potential  use  of  geotextiles  for  airfield  runways  is  considered.  The  conclusions 
of  Reference  13,  which  are  appropriate  to  the  present  research,  will  be  dis¬ 
cussed  in  some  detail  later  in  this  report. 

To  date,  most  of  the  theoretical  research  on  earth  reinforcing  has  been 
directed  towards  “classical"  reinforced  earth  retaining  walls  (References  14  and 
15).  Noteworthy  in  this  regard  is  the  work  described  in  References  16  and  17. 
Very  little  theoretical  research  has  been  carried  out  on  the  problem  of  rein¬ 
forced  embankments  or  pavements,  with  the  exception  of  the  work  described  in 
References  18,  19  and  20.  For  example, in  Reference  20, it  is  shown  that  horizon¬ 
tally  lying  reinforcement  layers  may  not  be  the  optimum  orientation  for  the 
reinforcement  under  embankments.  It  is  not  known  whether  a  similar  situation 
exists  with  respect  to  pavement  systems,  although  practical  construction 
requirements  probably  would  control  design. 

Related  theoretical  analyses  and  model  footing  tests  on  horizontally  rein¬ 
forced  foundation  soils  were  conducted  by  Binquet  and  Lee  (References  21  and 
22).  Their  results  indicated  that  the  settlement  behavior  and  ultimate  bearing 
capacity  were  improved  over  unreinforced  soils  by  factors  of  two  to  four.  The 
actual  improvement  depended  on  the  percentage  of  total  area  covered  by  the  rein¬ 
forcement  and  the  thickness  and  spacing  of  the  reinforcing  layers.  Increasing 
the  amount  of  reinforcing  definitely  improved  the  ultimate  bearing  capacity  and 
reduced  the  settlements,  especially  if  the  reinforcing  began  near  the  bottom  of 
the  Tooting.  This  result  is  similar  to  that  found  by  Haliburton,  et  al .  (Refer¬ 
ences  12  and  13).  Similar  results  were  also  obtained  in  some  of  the  present 
theoretical  analyses.  Improvement  in  bearing  capacity  resulted  even  when  the 
reinforcing  was  located  significantly  deeper  than  the  lowest  point  on  the 
theoretical  bearing  capacity  failure  surface,  an  observation  which  suggests  that 
a  different  failure  mechanism  is  operative  when  reinforcing  is  used.  This 
result  was  not  verified  previously  by  Haliburton,  et  al .  (References  12  and 
13). 


finally,  the  results  of  the  studies  of  Binquet  and  Lee  (References  21  and 
22)  indicate  that  the  greatest  advantage  of  reinforcing  foundations  was  for 
short-term  construction  involving  heavy  loads  over  inferior  foundation  condi¬ 
tions.  To  some  extent,  such  conditions  involve  the  function  of  separation 


rather  than  reinforcement.  However,  the  research  by  Binquet  and  Lee  was  limited 
to  smooth  noncontinuous  metallic  strips  (rather  than  continuous  sheets  of  other 
materials  as  tested  by  Hal i burton  and  in  the  present  study)  and  only  one  soil,  a 
sand,  was  tested.  Only  static  surface  loadings  were  applied  and  no  assessment 
was  made  of  the  possible  effects  of  impact  or  vibratory  loading  on  the  behavior 
of  the  reinforced  soil  system. 

2 .  RELATED  WORK  ON  REINFORCED  ROADS;  MESL 

Considerable  related  research  has  been  conducted  on  thin  but  relatively 
strong  tensi le- resistant  materials  or  membranes  under  small  embankments  (haul 
roads)  constructed  on  very  soft  foundations.  Much  of  this  research  has  been 
sponsored  by  the  manufacturers  of  nonwoven  geotextiles,  and  unfortunately  in 
most  cases,  the  assumptions,  theoretical  analyses,  and  backup  experimental  data 
are  not  publicly  available  (References  23  and  24).  The  procedures  developed  by 
Barenberg  and  his  students  (References  25,  26,  and  27),  which  are  based  on 
research  which  was  initially  privately  sponsored,  are  an  exception.  Kinney 
(1979)  has  developed  a  "fabric  tension  model"  by  which  the  modulus  of  the 
geotextile,  as  well  as  subgrade  strength,  traffic  loads,  and  rut  geometry  can  be 
appropriately  considered.  Summaries  of  the  various  manufacturer-sponsored 
design  methods  are  given  in  References  23,  24  and  28. 

Giroud  and  Noiray  (Reference  29)  developed  a  method  with  a  very  sound 
theoretical  basis  that  takes  into  account  full-scale  test  data  developed  at  the 
U.S.  Army  Engineer  (U.S.A.E.)  Waterways  Experiment  Station.  The  method  offers 
design  charts  that  allow  the  determination  of  aggregate  thicknesses  for  unpaved 
roads  when  geotextiles  are  used  as  reinforcement  and  when  traffic  is  taken  into 
account.  The  rut  depth  considered  in  the  design  charts  is  approximately  1  foot 
(0.30  meter).  Recently,  Sivakugan  (Reference  30)  has  prepared  design  charts  for 
lesser  rut  depths.  Important  for  the  present  research  is  the  fact  that  the  tire 
inflation  pressures  and  axle  loads  presented  in  the  charts  developed  by  Giroud 
and  Noiray  are  for  typical  construction  equipment  and  not  military  aircraft. 
Their  standard  axle  load  is  80  kN  (about  18,000  pounds)  and  maximum  tire  infla¬ 
tion  pressure  is  620  kPa  (about  90  psi).  Consequently,  considerable  extrapola¬ 
tion  would  be  necessary  to  use  their  charts  for  design  purposes  for  typical 
fighter  aircraft  loads  which  might  be  applied  to  alternate  launch  and  recovery 
surfaces.  Such  an  extrapolation  could  be  dangerous  without  the  results  of 
full-scale  tests  such  as  outlined  in  Section  V.  For  the  purpose  for  which  it 
was  developed,  the  Giroud  and  Noiray  method  is  simple  to  use,  and  according  to 
Giroud  (personal  communication,  1982)  the  method  has  been  used  with  considerable 
success  in  practice. 

Hamilton  and  Pearce  (Reference  31)  developed  guidelines  for  the  design  of 
flexible  pavements  using  slit  film  woven  fabrics.  The  method  is  specifically 
applicable  to  the  Texas  Gulf  Coast  region  where  very  poor  subsoils  predominate 
and  suitable  construction  aggregates  are  either  nonexistent  or  of  poor  quality. 
Significant  haulage  distances  can  result  in  extremely  high  construction  costs. 
Hamilton  and  Pearce  found  that  high  modulus  geotextiles  have  the  potential  of 
solving  many  of  the  pavement  problems  in  that  region.  They  present  a  design 
method  and  suggest  that  the  use  of  woven  geotextiles  offers:  (1)  a  reasonable 
and  cost-effective  alternative  to  mechanical  or  chemical  subgrade  stabilization, 
(2)  a  reduction  in  required  base  thicknesses,  and  (3)  an  extended  pavement  life. 


For  about  15  years,  the  U.S.A.E.  Waterways  Experiment  Station  (WES)  has 
been  conducting  research  on  membrane  encapsulated  soil  layers  (MESL).  Both 
nonwoven  geotextiles  as  well  as  military  membranes  such  as  the  T-16,  T- 17  and 
WX-18  membrane  mats  have  been  used  (References  32  and  33).  WES  has  also  con¬ 
ducted  research  on  "sand  bag"  type  structures  for  the  expedient  construction  of 
bridge  piers  and  abutments  in  a  theater  of  operations  (Reference  34).  Recent 
research  by  the  same  group  has  involved  very  unconventional  reinforcing  materi¬ 
als  such  as  small  plastic  cylinders  and  paper  and  aluminum  grids,  both  hexagonal 
and  rectangular  in  shape.  Roadways  were  constructed  of  reinforced  sand  on  both 
very  soft  clay  and  sand  subgrades.  These  tests  are  described  in  References  35, 
36  and  37 . 

3.  HAL  I BURTON 'S  WORK  AND  CRITIQUE  BY  VAN  DEN  BERG 

As  mentioned  above,  the  work  by  Haliburton,  et  al .  (References  12  and  13), 
is  particularly  pertinent  to  the  research  reported  herein.  This  section  will 
summarize  their  research  and  show  that,  in  some  cases,  their  conclusions  were 
verified,  and  in  other  cases,  they  were  not.  They  determined,  among  other 
things,  that  the  potential  improvement  of  the  performance  of  embankments  and 
runways  results  from  three  different  phenomena:  (1)  The  geotextile  appears  to 
act  as  a  separation  medium  which  prevents  the  intrusion  and  deterioration  of  the 
aggregate  materials  in  the  embankment.  This  phenomena  is  especially  pertinent 
when  the  subgrade  is  soft  and  cohesive.  (2)  There  appears  to  be  a  degree  of 
lateral  restraint  provided  by  the  fabric  to  the  cohesionless  cover  material.  As 
the  wheel  load  is  applied,  the  embankment  material  tries  to  spread  laterally  and 
this  expansion  tends  to  be  prevented  by  the  friction  developed  between  the 
fabric  and  the  embankment  material.  The  effect  is  to  increase  the  deformation 
modulus  of  the  material  in  the  embankment.  (3)  The  third  mechanism  or  benefit 
provided  by  the  geotextile  is  that  of  membrane-type  support.  For  this  mechanism 
to  occur,  relatively  large  deformations  must  occur  in  the  subgrade  to  mobilize 
the  full  membrane  resistance. 

Because  the  presence  of  the  fabric  is  the  controlling  factor,  there  is  an 
optimum  location  of  the  fabric  to  provide  for  lateral  restraint  reinforcement  in 
cohesionless  soils.  Haliburton  and  coworkers  suggest  that  this  location  is 
approximately  0.5  B  tan*  below  the  soil  surface,  where  B  is  the  effective  width 
of  the  loaded  area  and  *  is  the  angle  of  internal  friction  of  the  material  above 
the  fabric.  Placement  of  the  fabric  at  the  optimum  depth  produces  three  main 
types  of  behavior:  (1)  Elastic  deformation  of  the  cohesionless  soil  mass  above 
the  fabric  occurs,  with  considerably  greater  deformation  resistance  than  exists 
for  an  unreinforced  soil.  A  significantly  greater  resistance  to  initial  shear 
failure  is  also  developed.  (2)  After  the  initial  shear  strength  of  the  rein¬ 
forced  soil  is  exceeded,  rapid  yielding  occurs  and  the  loaded  area  sinks  into 
the  subsoil.  The  yielding  occurs  as  a  result  of  plastic  deformations  in  a  modi¬ 
fied  radial  shear  zone  around  the  loaded  area.  There  is  also  a  concurrent  shear 
failure  in  the  soil  mass  along  some  failure  surface  above  the  fabric.  (3)  The 
fabric  begins  to  "reinterfere"  with  the  plastic  flow  in  these  modified  radial 
shear  zones  as  sinkage  of  the  loaded  area  takes  place.  The  writers  believe  that 
these  zones  return  to  elastic  equilibrium,  which  stops  the  "sinking"  from  con¬ 
tinuing.  This  phenomenon,  a  “soil  strain  hardening,"  results  in  a  significant 
second-phase  strength  gain. 


Hall  burton  and  coworkers  found  that  if  the  fabric  is  placed  below  the 
optimum  depth,  only  minimal  improvement  in  the  load  deformation  characteri sties 
of  the  reinforced  soil  mass  will  be  realized.  However,  if  initial  soil  failure 
does  occur,  then  the  loaded  area  will  sink  to  near  the  optimum  depth  and  the 
strain-hardening  phenomena  will  begin  to  take  place.  In  practice,  such  sinkage 
would  likely  be  seen  as  excessive  rutting  in  the  wheel  paths  and  probably  would 
not  be  tolerable.  They  found  that  if  the  fabric  was  placed  at  the  optimum 
depth,  only  a  small  amount  of  loaded  area  sinkage  was  required  to  mobilize  a 
very  large  part  of  the  soil  strength  and  deformation  resistance  due  to  "strain 
hardening."  They  concluded  that  if  the  fabric  was  placed  on  a  "good  subyrade," 
then  the  second-stage  strength  would  likely  be  sufficient  to  carry  very  heavy 
imposed  loads,  although  no  data  were  presented  to  verify  this  belief.  If  a 
relatively  soft  subgrade  soil  is  present  under  the  fabric,  then  a  failure  in  the 
soil  under  the  fabric  is  likely  to  occur  which  would  prevent  full  utilization  of 
the  total  cover  soil  strength  available  from  the  fabric  reinforcement. 

In  contrast,  van  den  Berg  {personal  communication,  1981)  is  somewhat  criti¬ 
cal  of  the  laboratory  test  procedures  used  by  Haliburton  and  his  associates. 

His  primary  complaint  has  to  do  with  the  fact  that  the  fabric  was  clamped  at  the 
edge  to  a  frame;  that  is,  it  was  prestressed  to  some  extent  prior  to  burial  in 
sand  and  the  subsequent  loading  by  the  load  plate.  He  points  out  that  if  the 
failure  pattern  occurs  as  a  failure  wedge  in  the  classical  sense  as  postolateu 
by  Haliburton,  then  there  will  be  a  movement  outward  of  the  soil  mass  on  noth 
sides  of  the  failure  wedge.  The  resulting  shear  forces  will  tend  to  increase 
the  tension  forces  in  the  fabric.  Thus,  the  effects  of  "lateral  restraint"  and 
membrane  tension  will  be  difficult  to  separate.  Clamping  the  fabric  results  in 
a  fixity  that  is  almost  perfect,  and  van  den  Berg  is  not  sure  that  in  practice 
the  fabric  would  be  "anchored  in  that  perfect  way."  He  is  also  critical  of  the 
use  of  the  "low  quality"  base  materials  (sand)  and  the  fact  that  the  width  of 
the  fabric  is  only  six  to  eight  times  the  width  of  the  loading  plates.  He  con¬ 
cluded  that  there  may  be  something  like  an  "optimum  depth"  of  fabric  reinforce¬ 
ment  in  road  construction  but  the  concept  is  not  clearly  developed  nor  com¬ 
pletely  accepted  by  the  profession. 

4.  MATERIALS  FOR  REINFORCING 

Several  materials  have  been  suggested  as  being  potentially  useful  as  pri¬ 
mary  reinforcing  elements  in  airfield  runways  and  other  similar  construction. 
Excellent  descriptions  of  engineering  fabrics  or  geotextiles,  as  they  are  now 
called,  are  found  in  several  references,  particularly,  References  23  and  38. 

8oth  these  books  describe  geotextiles  in  some  detail  and  give  extensive  refer¬ 
ence  lists. 

New  materials  called  "geogrids"  have  recently  been  developed  which  have 
some  features  similar  to  geotextiles.  Geogrids  look  like  nets  of  plastic, 
although  the  strength  of  these  materials  is  significantly  greater  than  typical 
plastic  nets.  An  English  manufacturing  company,  Netlon,  produces  the  new 
material  under  the  trade  name  of  "Tensar®."  Tensar*'  geogrids  were  utilized  for 
a  few  of  the  tests  in  the  present  study,  and  their  properties  will  be  described 
in  some  detail  later  in  this  report.  On  a  weight  basis,  Tensar1  nets  are 
stronger  than  steel  but  their  costs  are  on  the  order  of  woven  geotextile  materi¬ 
als.  They  also  havo  an  added  advantage  over  ordinary  geotextiles  of  providing 


"interlock,"  in  addition  to  frictional  resistance,  if  materials  coarser  than 
sands,  e.g.,  fine  gravel,  are  used  in  construction. 

5.  RELATED  WORK  ON  RRR  AND  ALRS 

Much  of  the  research  on  rapid  runway  repair  procedures  has  been  summarized 
in  References  39  and  40.  In  Reference  41,  materials  for  "contingency  runways" 
were  evaluated  by  full-scale  tests  with  the  F-4  load  cart.  An  open  graded 
crushed  stone  base  course  with  a  sod  surface  was  found  to  rut  severely  when  sub¬ 
jected  to  the  F-4  loadings.  Some  success  has  been  reported  with  cement  and  lime 
stabilization,  but  almost  any  kind  of  unstabilized  soil  or  sod  surface  (overly¬ 
ing  the  stabilized  base)  was  found  to  rut  excessively  unless  the  top  surface  was 
less  than  1-inch  thick.  With  such  thin  surface  layers,  the  topsoil  might  be 
subject  to  sliding  on  the  base  course  due  to  the  breaking  action  of  aircraft 
tires.  It  is  possible  that  geotextiles  might  improve  the  bonding  of  a  thin  top¬ 
soil  surface  to  stabilized  subbases.  Furthermore,  we  suggested  that  the  top 
layer  of  contingency  runways  and  alternate  launch  and  recovery  surfaces  might  be 
stabilized  economically  and  reinforced  by  geotexti le-1 ike  materials.  The  use  of 
geotextiles  on  ALRS  sites  in  conjunction  with  vegetation  established  on  the  sur¬ 
face  would  enhance  their  multiple  use  and  add  considerably  to  the  concealment  of 
the  site. 

6.  OBJECTIVES  OF  THE  PRESENT  RESEARCH 

A  major  objective  of  the  present  research  was  to  develop  an  improved 
theoretical  model  for  the  interaction  of  reinforcement  in  soil.  Such  a  model 
would  improve  our  understanding  of  the  fundamental  interaction  of  subgrade  soils 
with  various  types  of  reinforcement  and  pavement  materials  when  subjected  to 
surface  loadings.  An  additional  objective  was  to  develop  reinforcement  confi¬ 
gurations  for  ALRS  pavement  systems  and,  if  possible,  optimize  such  systems. 
Finally,  plans  were  to  be  developed  for  full-scale  field  tests  to  be  carried  out 
as  "Phase  II." 

7.  OUTLINE  OF  THIS  REPORT 

In  Section  II  of  this  report,  the  analytical  work  carried  out  will  be  sum¬ 
marized.  In  Section  III,  the  details  of  the  experimental  phase  of  the  research 
will  be  given  and  the  results  of  the  experiments  will  be  examined  in  terms  of 
the  theoretical  developments  and  empirical  physical  behavior  models.  A  brief 
section  is  also  included  on  the  implications  of  the  present  research  for  RRR  and 
ALRS.  Recommendations  for  Phase  11  research,  field  studies,  are  also  outlined. 
Finally,  some  conclusions  and  recommendations  are  presented. 


SECTION  II 


ANALYTICAL  WORK 


1.  INTRODUCTION 

It  is  the  function  of  analysis  in  Civil  Engineering  to  perform  inexpensive 
experiments  (generally  on  paper)  that  will  help  to  predict  the  response  of  other¬ 
wise  expensive  prototype  systems.  Considering  the  number  of  possible  permuta¬ 
tions  of  soils,  types  and  placement  of  fabrics,  pavement  types,  and  aircraft 
gear  configurations  and  loading  involved  in  this  study,  it  is  apparent  that 
discovered  truths  must  first  be  approached  by  means  of  analysis.  On  the  other 
nana,  the  dominating  factors  in  the  analysis  of  an  aircraft  moving  over  a  pave¬ 
ment  system  are  extremely  complex  and  are  generally  "wished  away"  or  ignored  in 
ordinary  mechanics.  Consider  this:  the  pavement  system  presents  itself  as  a 
number  of  layers;  each  one  is  thin  in  comparison  to  the  contact  dimensions  of 
the  load  (the  tire  imprint).  The  layers,  in  turn,  are  composed  of  discrete  par¬ 
ticles  that  are,  at  best,  cemented  by  a  thermal-sensiti ve  material.  Tne  system 
may  exhibit  varying  degrees  of  saturation  by  liquids  and/or  gases.  Loadings  are 
xnown  only  in  magnitude,  not  in  their  distribution.  Contacts  between  layers  vary 
with  time,  ambient  conditions  and  loadings. 

In  spite  of  the  noted  complexities,  the  engineering  system  offers  a  method¬ 
ology  witn  a  high  potential  for  success;  namely,  (a)  to  offer  simple  analytical 
models  that  relate  the  pertinent  factors  believed  to  govern  performance  and  (b) 
to  conduct  simple  laboratory-scale  model  tests  to  examine  the  various  offered 
analytical  developments. 

Analytical  and  experimental  models  may  themselves  be'  divided  into  three 
parts:  (a)  input,  (b)  transfer  mechanism,  and  (.;')  output.  If  both  types  of 
models  are  subject  to  the  same  input,  the  test  of  the  adequacy  of  an  analytical 
me  del  will  oe  its  ability  to  predict  the  observed  response  of  its  experime^ta1 
counterpart.  Several  such  models  will  be  presented  in  this  chapter.  All  models 
consider  the  loading  to  act  normal  to  the  plane  of  the  surface. 

7. .  HOMOGENEOUS  MEDIUM 

a.  plane  Strai  n  Case 

This  solution  was  presented  oy  narr  (Reference  42)  .  The  expected 
value  of  vertical  stress  (S  )  under  a  line  load  of  intensity  r,  per  unit  run, 
(Figure  I)  is  given  by  1 


i .  TTarr  (Reference  4?,  Chap.  7!  demonstrated  that  the  stress  at  a  point  in  a 
particulate  media  is  a  Pcisson  variable.  Consider  that  at  such  a  point 
there  may  exist  -  a  void  -  a  solid  -  a  liquid  -  a  gas  -  ...  -  and,  it  is 
evident,  that  any  measure  of  the  distribution  of  boundary  energy,  such  as 
stress,  at  a  point  is  indeed  a  variance.  The  expected  value,  the  meat 
value,  is  the  best  measure  of  this  intensity  at  a  point. 


1 


Figure 


1.  toad  Geometry. 


Harr  called  v  the  "coefficient  of  lateral  earth  stress."  He  showed  it  to  be 
related  to  the  more  common  "coefficients  of  lateral  earth  pressure."  The  clas¬ 
sical,  two-dimensional,  elastic  solution  (Flamant  solution)  requires  no  measure 
of  the  properties  of  the  material.  The  attenuation  of  stress  is  assumed  to  be 
purely  geometric. 

For  a  uniform  normal  load,  q,  over  an  infinite  strip,  (Figure  2),  of  width 
a,  the  equivalent  solution  is 


z» 

where  )  is  the  well -tabul ated  cumulative  distribution  (see  Reference  42,  p. 

459) . 

b.  Three-Dimensional  Case 

For  a  single  concentrated  force,  P,  normal  to  the  plane  of  the  surface, 
Reference  42  gave  for  the  expected  value  of  vertical  stress 


This  expression  is  in  the  same  form  of  that  given  by  Boussinesq.  However, 
the  above  acknowledges  the  characteristics  of  the  media  (v);  whereas,  the  Bous¬ 
sinesq  solution  is  devoid  of  any  material  properties  characteristics. 


For  uniform  normal  load,  q  .over  a  circular  area  (Figure  3),  Equation  (3) 
generalized  to 


xr  i  or.  ■'•)  't  dr-  d 


general  solutions  of  equation  (4)  can  only  be  obtained  by  numerical  means, 
bone  results  will  be  presented  below. 

For  the  special  case,  under  the  center  of  a  uniformly  distributed  load.q, 
acting  on  a  circular  area  of  radius  a,  the  expression  reduces  to 


;.  multilayer  MEDIUM 

Altnouyh  many  solutions  exist  in  the  literature  for  the  transmission  of 
stress  through  a  layered  system,  none  are  unique,  even  with  the  constraints  of 
the  classical  theory  of  elasticity.  The  calculated  values  of  stress  depend  on 


Figure  2.  Infinite  Strip  Geometry. 


Figure  3.  Uniform  Load  over  Circular  Area. 


the  assumptions  invoked  at  the  boundaries  between  the  individual  layers.  All 
solutions  assume  continuity  of  di spl acements  across  interfaces;  however,  this 
assumption  alone  is  not  sufficient.  Some  statement  must  be  offered  as  to  the 
transformation  of  shearing  stresses  at  the  interfaces  as  this  is  not  known  a 
priori.  Two  such  assumptions,  commonly  prescribed,  are  that  (a)  there  is  "no 
slipr"  and  ib)  there  is  “no  friction''  at  the  interface.  Although, ii  ir 
era! ly  believed  that  these  provide  "upper"  and  “lower"  bounds  to  the  system; 
unfortunately,  this  is  not  the  case.  A  boundary  wherein  one-half  evidences  no 
sMp"  and  the  other  half  the  "no  friction"  condition  would  be  more  severe  than 
either.  Harr  (Reference  4c)  presented  a  solution  that  obviates  the  need  to 
specify  the  shearing  stresses  at  the  interfaces  of  layered  media.  For  the 
expected  vertical  stress  due  to  a  line  load  of  intensity  p  (Figure  4),  he  gave 
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is  the  vertical  distance  into  the  Nth  layer. 
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As  an  trample;  31. EN  A  three-layer  system  is  subjected  to  a  line  load 
ot  ijCJ  lb  ft  with  the  id  lowing  information:  hj  -  i  ft,  v,  =  0.4;  -  i  ft, 

.  =  u.j.  ;  :j.4,  tv  -  unbounded.  Find  the  expected  value  for  the  vertical 
nbrnal  stress  a  feet  into  the  third  layer  immediately  under  the  line  load. 


SOLUTION  From  equation  (6b),  the  equivalent  thickness  is 


>us,  from  Equation  (ba)  we  have  for  the  expected  vortical  normal  stress  at  a 
depth  of  3  feet  in  the  thirq  layer,  immediately  under  the  line  load  (x  -  o'. 


Ine  theory  of  elasticity  (Reference,^)  would  give  for  this  case  (assuming 
a  homogeneous  section),  -  HS4.9  lb,  ft  . 

^or  a  uniformly  distributed  normal  load  q,  acting  over  a  circular  area  of 
’■dd’us  a,  the  equivalent  form  of  Equation  16)  is 


Figure  A 


Again,  results  can  only  be  obtained  by  numerical  methods,  examples  will  be 
developed  subsequently. 

4.  SOIL -FABRIC  MEDIUM 

An  anaytical  model  was  developed  tc  simulate  the  action  of  a  soil -fabric 
system.  Briefly,  the  model  assumed  a  particulate  medium,  as  described  i>y  equa¬ 
tions  (I)  through  (?},  founded  on  a  layer  specified  as  a  Winkler  body.  Between 
the  two  layers  lat  their  interface)  is  a  "geotextile"  membrane.  Solutions  were 
obtained  oy  numerical  means  for  a  range  of  values  of  soil  and  yeotextile  proper¬ 
ties.  The  details  of  the  model  and  some  results  were  given  in  a  paper  by  the 
authors  and  Mr.  P.  L.  bourdeau,  which  was  presented  at  the  Second  International 
conference  on  ueotext i les ,  August  1982  (Reference  44).  A  copy  of  the  paper  is 
in  Appendix  A.  The  results  of  the  developed  analysis  indicated  that  very  little 
improvement  in  the  load-carrying  capacity  of  the  reinforced  system  would  be 
achieved  unless  the  subyrade  was  very  soft.  f:his  can  be  seen  in  Figures  w  and 
lo  in  the  paper  in  Appendix  A.!  Another  interesting  observation  was  that  the 
effective  length  of  the  reinforcement  was  relatively  snort.  In  other  words,  at 
a  rather  short  distance  from  the  loaded  area,  the  re i  r. forcenent  would  not  'feel'1 
any  stressing  due  to  the  surface  load.  This  result  seemed  anomalous  at  first, 
but  recent  evidence  by  Andrawes,  et  al.  (Reference  45),  showed  very  similar 
results  for  plain  strain  tests  on  reinforced  sands.  They  found  that  the  meas¬ 
ured  tensile  force  in  the  rei of orcement  was  practically  nonexistent  at  1:  tances 
from  the  centerline  as  close  as  twice  the  diameter  of  width  of  the  loade-1  area. 

Although  the  i.ioael  does  not  consider  failure  induced  by  insufficient  bear¬ 
ing  capacity,  it  does  investigate  two  other  important  possible  modes  of  dis- 
t re. a:  u)  that  of  the  tensile  force  exceeding  the  tensile  strength  of  the 
fabric,  and  (?)  the  slipping  ("pull  out")  of  the  fabric  if  its  length  is  less 
than  that  required  for  stability. 

-  .  STRAIS  -JifFuilON 

the  previous  cases,  as  in  continuum  mechanics,  presume  that  the  stress  dis¬ 
tribution  at  the  surface,  the  input  into  the  system,  is  a  known  quantity.  This 
■ >  seldom  the  case,  contact  pressures  under  tines  are  tar  from  the  idealised 
uni  term  values  commonly  assumed,  '.n  addition,  the  pressures  change  with  time, 
temperature  ana  ambient,  conditions  (Reference  4b). 

The  rationale  of  continuum  mechanics  is  to  determine  the  transmission  ot 
boundary  energy  through  media  as  statements  of  stress  intensities.  Granted  such 
measures,  constitutive  relationships  are  then  introduced  that  relate  strides  to 
commensurate  strains.  Th.e  vertical  components  of  the  strains  are  then 
integrated  to  obtain  the  surface  deflections,  which  are  examined  relative  to 
differential  settlements  and/or  rutting. 

When  one  considers  the  mechanism  responsible  for  t.he  deformations  noted  at 
the  su  rface  of  a  medium  subject  to  induced  loading,  it  is  apparent  that  it  is 
not  caused  oy  the  deformation  of  the  individual  particles,  tven  a  cursory  exam¬ 
ination  will  demonstrate  that  the  registered  displacements  occur  in  response  to 
the  relative  movements  of  the  more  mobile  individual  particles  to  form  a  denser 
matrix.  Motivating  this  movement  is  the  induced  boundary  energy.  in  principle. 


when  a  rigid  plate  is  impressed  upon  a  soil  surface,  the  medium,  in  effect, 
"diffuses"  this  displacement. 

The  question  arises,  if  boundary  strains  (or  displacements)  are  known  pre¬ 
cisely  and  if  we  seek  their  diffusion  into  media,  why  not  do  this  directly? 

This  suggests  a  "strain  diffusion"  process  rather  than  the  more  conventional 
stress  diffusion  process.  The  investigation  of  this  matter  was  considered  to  be 
germane  to  the  present  study  and  hence  was  conducted  as  an  integral  part  of  the 
work  project.  Its  development  follows. 

Adopting  the  foregoing  diffusion  of  stress  as  a  strain  process,  as  was  done 
to  obtain  Equation  (7),  the  expected  vertical  strain  at  any  point,  P,  due  to  a 
surface  strain  of  intensity  q,  acting  over  a  circular  area  of  radius  a  (Figure 
is 

3  '  i  P  ^ 

—  ,  q  x‘  +  r”  -  ?xr  cos  r>  , 

?-vz‘  o  0 

where  u  is  a  material  parameter.  Defining  p(x)  as  the  deflection  of  the  sur¬ 
face,  we  have 

cf' 

o(x!  -  /  (x,z)  d?. 

r\  •  • 

introducing  the  Bessel  functions  of  the  first  kind,  J^(x)  and  J  (x),  Equa¬ 
tion  (b)  takes  the  form 

no  ..2 

•  _  -  r.  j  J  'yt!  J .  (at)  exp  i-  i’ 

where : 

•n 

t  in 

j  .  *  ;  :  '  4 ,  ‘  •  ;r 
rl  '  kO 

r or  1 arge  val ues  of  t, 


Equations  (9)  and  (10)  were  solved  numerically.  In  Equation  (lu),  t  is  a 
dummy  variable  introduced  by  the  transformation  t  =  t(x,  a).  It  should  be  evi¬ 
dent  tnat  Equation  (10)  lends  itself  to  simpler  numerical  solutions  than  dees 
Equation  (8).  Examples  will  be  given  below. 

Finally,  a  parallel  situation  to  that  in  Section  II. 4  on  "Soil-Fabric 
Medium"  was  investigated  using  "strai n-diffusion . "  Equations  (9)  and  (10)  were 


assumed  to  hold  for  each  of  the  layers  (Figure  t>) .  For  the  multilayar  case, 
that  is  of  (h.,u,l  and  (h„,^),  Iquation  (10)  must  be  solved  twice  for  each 
value  of  the  parameter  u,  i.e.. 


tr:d  >  x , 1  I 

Iquation  (9)  then  takes  the  form  of 


o .  EXAMPLES  OF  MODELS 

a.  Homogeneous  Medium  -  Infinite  Strip 

Figure  7a  shows  a  plot  of  the  expected  vertical  stress  for  a  unit  Jis- 
triduted  load  on  an  infinite  strip.  Figure  2.  Equation  (2)  governs  this  case. 
The  depth  is  taken  as  one-half  the  width  of  the  strip  (z  =  a),  the  "v  parameter 
is  0.2S,  and  a  unit  load  is  applied. 

b.  Homogeneous  Medium  -  Circular  Area 

Figure  7b  provides  a  plot  of  the  vertical  stress  under  the  center  of  a 
uniformly  loaded  area  of  radius  a  at  a  depth  equal  to  the  radius.  Equation  !v' 
applies  for  this  case.  The  differences  between  Figures  7a  and  b  are  only  minor. 
That  is,  only  a  small  part  of  the  energy  available  for  the  infinite  strip  is 
effective  under  the  center  (or  centerline)  of  the  loaded  area. 

<: .  Strai  n  Diffusion 

Recall  in  Section  II. 5  that  a  methodology  was  developed  that  invokes 
the  diffusion  of  boundary  di spl acements  through  a  continuum  rather  than  throug1 
the  customary  boundary  stresses ■  ThTierent  in  the  new  analysis  is  the  u- 
parameter.  Equation  (A),  which  scales  the  transmission  of  the  boundary  energy. 

A  number  of  numerical  values  were  assumed  for  this  parameter  to  examine  the 
influence  that  it  had  on  surface  deflections.  It  was  assumed  that  a  uniform 
unit  displacement  acted  over  a  circular  area  of  radius  a;  that  is.  Equation  i-' 
applied  with  q  =  1.  The  relative  depth  of  the  medium  was  h  =  7./  (the  actual 
depth  in  the  test  box  of  2d  inches  divided  i>y  a  plate  diameter  of  3.0  inches). 

A  computer  program  (Program  STRAIN)  for  calculating  the  surface  deflections  was 
written  and  is  given  in  Appendix  B.  Input  to  this  program  were  the  values  of  e 

=  1,  ar  1,  h-7.7,  and  various  values  of  u,  ranging  from  3.0  to  nearly  that  of 

zero  i D.UOS).  Output  was  in  the  form  of  normalized  tabulated  values  of  deflec¬ 
tion  (actually  a  deflection  ratio  equal  to  the  deflection  at  any  radius  divided 
by  the  maximum  deflection)  and  the  relative  radius,  x/a.  The  results  for  a 
number  of  values  of  m  are  plotted  in  Figures  8  to  14.  As  an  example  of  how 
these  r lots  were  obtained,  the  output  for  cases  m  -  0.3  and  t.  =  0.2  are  tabu¬ 
lated  :,low  and  plotted  in  Figures  11  and  12,  respectively. 
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Figure  8.  Relative  Surface  Deflections  Given  a  Uniform 
Displacement  over  3  Circular  Area  of  Radius 
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Figure  11.  Relative  Surface  Def 
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gure  1 A ,  Relative  Surface  Deflections  Given  a  Uniform 
Displacement  over  a  Circular  Area  of  Radius 
A,  4  -  0. 
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it  was  desired  to  obtain  the  value  of  the  parameter  that  agreed  with 
experimental  observations  of  the  deflected  shape  of  the  surface.  Since  this 
required  a  "back-calcul  ation:'  of  u  for  the  soil,  no  fabric  was  assumed  to  exist 
for  this  case.  It  is  seen  from  figures  8  to  14  that  as  u  decreases  (figure 
numbers  increasing),  more  heave  is  evident.  The  heaving  of  the  soil  is  of  con¬ 
siderable  importance,  as  it  was  observea  in  all  the  laboratory  tests  (see  Fig¬ 
ures  0 1  and  62)  to  some  degree.  It  should  be  emphasized,  as  shown  in  figure  '. 
t»e  das-ical  theory  of  elasticity  cannot  account  for  the  heaving  of  a  nomocene- 
ous  nedium. 


Figures  61  and  o2  present  normalized  deflection  basin  measurements  for  a 
number  of  physical  tests.  It  is  seen  that  heaving  on  the  order  of  20  percent  of 
the  induced  deflection  is  not  out  of  order.  Comparison  with  Figures  i>  to  14 
indicates  that  the  ..-parameter  is  very  close  to  zero.  However,  „  =  u  introduces 
mathematical  instability  in  liquation  (8).  Consequently,  a  value  of  G.Gu5  was 
used  vi  subsequent  work.  The  physical  significance  of  this  parameter  is  diffi¬ 
cult  tc  ascertain  at  this  time.  Apparently,  it  is  somewhat  analogous  to 
Poisson's  ratio. 

Once  the  ..  parameter  for  uniform  soil  (^ }  is  computed,  it  can  also  be  used 
to  obtain  the  „  parameter  fc->  foe  fabric  (i.e!,  Some  calculations  were 

made  to  determine  the  value  of  ...  using  the  experimental  results  for  the  cirr  'ar 
plate  (:igjre  61).  Then  the  calculated  values  of  ^  and  were  used  in  Equa¬ 
tion-  > 9)  and  ,10)  to  predict  the  surface  deflection  pattern  of  a  soil-fabric 
syste"'  when  loaded  in  plane  strain.  "The  predictions  did  not  agree  well  with  the 
experimental  results. 


d.  Displacements  and  Settlements 

The  objective  of  all  the  foregoing  analyses  is  to  predict  the  deflected 
shaoe  of  a  surface  when  subjected  to  an  imposed  load.  After  this  objective  is 
achieved,  determinations  can  be  made  as  to  whether  the  degree  and  natu1'-  the 

defornations  can  be  tolerated  without  impairing  the  intended  use  of  the  facil- 
i  ty. 


As  noted  in  Section  b  and  Figure  7,  classical  mechanics  predict  stress 
distributions.  To  translate  these  into  surface  displacements  recourse  must  do 
nan  to  constitutive  relationships  between  stress  and  strain. 

As  is  commonly  the  case,  Hooke's  law  will  be  generalized  to 


where  "S  and  o  are  normal  stresses  on  planes  parallel  to  the  unloaded  surface 
and  £  afid  ,.  ar^  Young's  modulus  and  Poisson's  ratio,  respectively  (sec  Reference 
43  for  the  derivation!.  For  axi symmetri c  conditions.  Equation  (11a)  reduces  to 


where  3  and  S  are  the  vertical  and  radial  normal  stresses,  respective1.-, 
z  r 

Harr  fReferer.ee  4?)  found  that  the  radial  normal  stress  can  be  related  to 
the  vertical  normal  stress  by 


Substituting  this  expression  with  S  from  Equation  (5)  into  Equation  ; \ \ n ) 
produces  an  equation  for  the  strain  unde?  the  center  of  a  circular  load  with  a 
uni  form  normal  intensity  q,  or 


Equation  (ip)  can  tie  expressed  in  the  form 


"'•ere  I  (  ',  z  .  y)  is  a_  "vertical  strain  influence  factor"  that  is  a  fun. 
of  the  parameters  -  and  ...  and  z  is  the  depth/radius  ratio  of  a  circular  !  a 
ar“a  , z  -  z  a. 

r  quit i on  (dj  gives  the  displacement  of  the  surface  at  the  center  ot  tre 
ci'-cular  lead,  which  can  be  combined  with  Equation  U3)  to  give 


The  stress  approach  could  be  extended  to  the  plane  strain  condition  by 
merely  selecting  the  I  diagram  under  such  conditions  corresponding  to  the  com¬ 
puted  value  for  sand.  2 Modulus  plays  an  important  role  in  the  settlement  oorpu tui¬ 
tion,  and  because  of  the  dependency  of  modulus  on  the  degree  of  prestressing,  ar. 
experimentally  obtained  modulus  is  required. 
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SECTION  III 


EXPERIMENTAL  WORK 


1.  INTRODUCTION 

The  objectives  of  the  experimental  portion  of  the  subject  research  w.-ro  to: 

a.  Obtain  parametric  data  for  the  theoretical 
formulation  developed  in  Section  II. 

b.  Test  the  theoretical  formulation  developed. 

c.  Observe  the  empirical  behavior  of  fabric- 
reinforced  sands  to  obtain  an  indication  of  the 
relative  improvements  in  system  performance  due 
to  the  presence  of  the  geotextile. 

d.  Verify,  insofar  as  possible,  previous  research 
findings. 

2.  DESCRIPTION  Of'  EXPERIMENTAL  PROGRAM 

A  series  of  laboratory-scale  model  loading  tests  of  reinforced  and  unrein¬ 
forced  sanas  was  conducted  and  analyzed.  A  detailed  description  of  this  portion 
of  the  research  follows. 

a.  Variables 

As  suggested  previously,  the  number  of  possible  permutations  of  the 
properties  of  soil,  rei nf orcement ,  pavement  types,  and  aircraft  gear  configura¬ 
tions  and  loadings  is  very  large.  Thus,  to  maintain  the  size  of  the  experimen¬ 
tal  program  within  time  and  budget  constraints,  it  was  decided  to  hold  constant 
a  number  of  possible  variables  throughout  the  experimental  program.  For  exam¬ 
ple,  to  minimize  the  effect  of  soil  variability  as  much  as  possible,  only  one 
soil,  a  relatively  uniform  Ottawa  sand,  was  tested  at  a  constant  relative  den- 
s  i  ty , 


Sand  w3S  chosen  for  Phase  I  tests  because  it  is  easier  than  cohesive 
materials  to  handle  in  large  volumes  in  the  laboratory  and  to  control  densities, 
furthermore,  the  beneficial  effects  of  reinforcing  sands  is  known  from  previous 
research,  whereas  reinforcement  of  cohesive  soils  with  geotextiles  has  not  often 
beer.  done.  A  coarser  material,  e.g.,  a  fine  gravel,  could  have  been  used  to 
model  the  granular  surface  materials,  but  it  was  felt  that  a  better  indication 
of  the  degree  of  improvement  due  to  geotextile  reinforcement  would  be  shown  with 
a  poor  quality  surface  material  such  as  uniform  sand.  This  was  thought  to  pro¬ 
vide  a  "lower  bound"  of  granular  soil  performance.  A  single  relative  density 
(medium  dense)  was  chosen  as  being  typical  of  field  deisities  obtainable  with 
ordinary  construction  equipment. 

Ideally,  tests  should  have  been  conducted  using  elliptically  shaped  plates 
approximately  9  inches  in  average  diameter,  but  the  size  of  the  test  box,  3c 
inches  square,  was  a  constraining  •  or.  Since  the  scale  ratios  were  only  i.f 


and  3,  it  was  believed  that  the  results  of  tests  on  3-  and  b-inch  diameter  cir¬ 
cular  plates  could  be  extrapolated  to  provide  reasonable  estimates  of  the  per¬ 
formance  of  a  9-inch  diameter  loaded  area.  Rigid  plates  were  chosen  not  only 
for  convenience,  but  with  the  knowledge  that,  as  far  as  the  soil  was  concerned, 
with  aircraft  tire  pressures  of  the  order  of  250  to  300  psi ,  tire  contact  areas 
are  effectively  rigid.  In  addition  to  circular  plates,  a  few  tests  were  con¬ 
ducted  with  a  rigid  plate  in  a  plane  strain  configuration.  A  loading  plate  i 
indies  wide  fay  30  inches  long  (almost  the  entire  width  of  the  test  box)  was  used. 
These  tests  were  an  attempt  to  provide  experimental  data  to  verify  the  plane 
strain  theory  described  in  Section  II. 

for  the  first  series  of  tests,  a  woven  polyester  high-modulus  fabric  was 
chosen,  primarily  because  considerable  testing  had  already  been  carried  out  0" 
this  fabric  at  Purdue  University  and  elsewhere.  The  second  type  of  reinforce¬ 
ment  tested  was  high-strength  extruded  plastic  (polypropylene  and  polyethely ne ! 
geogrid  nets.  Geometric  variables  included  the  number  of  layers  of  reinforce¬ 
ment  and  the  depth  and  spacing  of  those  layers,  as  illustrated  in  Figure  15.  Tc 
provide  a  "common  denominator"  against  which  to  measure  the  degree  of  improve¬ 
ment  of  the  rei nforcement ,  several  tests  were  run  without  any  reinforcement. 

Table  1  lists  the  test  variables,  the  quantities  assigned  to  those  vari¬ 
ables,  and  the  symbol  adopted  for  each  variable. 

c.  Code 

To  facilitate  identification  of  the  individual  tests,  a  coding  system 
war<  adopted,  usi  ng  the  symbols  for  the  test  variables  listed  in  Table  1.  For 


ci rcul ar 

plate 

load  tests 

,  the  code  is 

number 

type  of 

edge 

( di ameter 

)(CP) 

of  layers 

reinf.  -  (depth)  - 

(spacing)  - 

cond 

0 

-- 

3 

CP 

1 

1 

W  d ,  i  n . 

0 

PF 

b 

CP 

2 

0 

s,in. 

FF 

If  an  item  is  omitted  for  one  reason  or  another,  the  symbol  is  left  blank.  For 
example,  the  symbols  for  the  tests  without  rei nforcement  are  3CP  and  6CP.  Those 
with  only  one  layer  of  reinforcement  omit  the  symbol  for  s,  and  so  forth. 

For  the  plane  strain  tests  without  fabric,  the  test  symbol  is  simply  PS. 

For  example,  with  one  layer  of  geogrid  reinforcement  at  2  inches  depth,  the  test 
designation  is  PS1G-2. 

Table  2  lists  the  tests  performed  and  their  symbols. 

3.  DESCRIPTION  OF  APPARATUS  ANU  EQUIPMENT 
a.  Loading  box  and  Reaction  Frame 

A  test  box  of  plywood  and  steel  angle  sections  was  constructed  for  this 
project.  The  inside  dimensions  of  the  box  are  80  centimeters  by  80  centimeters 
and  about  70  centimeters  high  (31.5  x  31.5  x  27  inches).  Because  of  clearance 


™BLE  1.  TEST  VARIABLES  AND  SYMBOLS 


V  a  r  1  abl  e 

1.  Soil  Type 

2.  Relative  Density 

3.  Loaded  Area 
Shape  and  size 


4.  Reinforcing 


I)  i  sp_0£i_tion_ 

Ottawa  "flint"  sand 
Medium  dense  (D  7(Ji) 


Circular  plate: 

3  in.  diameter 
6  in.  diameter 
Plane  strain: 

3  in.  x  30  in. 

Rigid 

Quasi-static  (approx, 
constant  rate  of  load) 


None 

Woven  polyester 
Oeogri ds 


Symbol 


Number  of  Layers, n 


Deptn  of  rein¬ 
forcing.  d 

spacing,  s 


Vari able 
Vari able 


(Depth  in 
inches) 

(Spacing  in 
inches) 


edge  conditions  r ree 

J  Partial  fixi 


Partial  fixity 
Full  fixity 


r 


TABLE  2. 

TESTS  CONDUCTED 

AMO  THEIR  SYMBOLS 

Total  tests  Tests  with 

Code 

Explanation 

Conducted 

Usable  Data 

3CP 

without  fabric 

2 

2 

3CP1W-G .5 

with  one  layer 

2 

2 

3CP1W-1.2  (or  1) 
3CP1W-2 

of  fabric 

2 

1 

1 

1 

3CP1G-1 

3CP1G-2 

with  one  layer 
of  SS2  grids 

1 

1 

3CP1W-1-PF 

with  partial 
edge  fixity 

2 

1 

3CP1W-1-FF 

with  full 
edge  fixity 

1 

1 

6CP 

without  fabric 

8 

4 

6CP1W-1 

with  one 

2 

2 

6CP1W-2 

layer  of 

2 

2 

6CP1W-3 

f abri c 

2 

2 

6CP2W-U  .8-0  .b 

with  two  layers 

1 

1 

1 

OLP2W-1 .6-1.6 

of  fabric 

1 

PS 

without  fabric 

1 

2 

PS1W-1 

with  one 

1 

1 

PS1W-2 

layer  of  fabric 

1 

1 

PS1G-1 

with  one  layer 

1 

1 

PS1G-Z 

of  SR2  grids 

1 

I 

j 


requirements  for  the  load  cell  and  the  loading  plate,  the  depth  of  sand  tested 
was  always  about  45  or  50  centimeters  (18  to  20  inches).  Two  0.5-inch  thick 
plywood  sheets  formed  the  sides  of  the  box,  which  was  reinforced  by  steel  angle 
sections  held  together  with  long  threaded  tie  rods.  Figure  16  is  a  photograph 
of  the  test  box. 

A  heavy  steel  test  frame  was  available  in  our  laboratories  from  previous 
research  on  model  asphalt  pavements  (Reference  47 ) .  The  reaction  frame  sup¬ 
ported  the  hydraulic  load  actuator  and  other  parts  of  the  loading  system,  as 
shown  at  the  top  of  Figure  16.  Figure  17  is  a  photograph  of  the  test  frame. 

b.  MTS  "Closed  Loop"  Hydraulic  Loading  System 

The  MTS  loading  system  used  in  this  research  was  manufactured  by 
Research,  Inc.,  of  Minneapolis,  Minnesota.  Figure  18  is  a  photograph  of  the 
controlling  system  (in  background):  the  loading  actuator  and  accumulator  were 
shown  on  the  top  of  the  load  frame  in  Figure  16.  The  entire  MTS  loading  system 
consists  of  the  following  components: 

Servoram  Hydraulic  Actuator,  Model  204.13,  5  kip, 

b-inch  stroke 

Servoram  Hydraulic  Actuator,  Model  202.03,  30  kip, 

b-inch  stroke 

Hydraulic  Power  Supply  (3000  psi).  Model  502.03 

Electronic  Control  System: 

Function  Generator,  Model  410.21 
Counter  Panel,  Model  417.01 
Control  Panel,  Model  413.04 
"Servac"  Control,  Model  401.02 
Recorder  Input  Selector,  Model  414.02 
Transducer  Conditioner  Panel,  Model  425.41 

The  MTS  loading  system  is  a  "closed  loop”  system  because  transducers  and 
electronics  in  the  system  control  the  pressures  and  deformations,  so  that, 
within  system  capability,  the  desired  load- deformation  time  history  can  be 
applied  to  the  specimen.  The  system  can  be  used  in  a  "stroke"  control  mode, 
wherein  the  deformation  of  the  actuator  piston  is  controlled  electronically;  or 
in  a  "load"  control  mode,  wherein  the  rate  and/or  magnitude  of  load  are  con¬ 
trolled  el ectroni cal ly .  The  stroke  mode  was  used  for  the  test  setup  when  "zero¬ 
ing"  the  load  plate,  etc.;  the  load  node  was  used  when  applying  the  load  to  the 
test  specimen. 

c.  Sand 

The  test  soil  selected  for  study  was  an  Ottawa  "flint"  silica  sand 
obtained  from  the  Ottawa  i  1  i c a  Co.  of  Ottawa,  Illinois.  This  sand  had  been 
used  in  previous  research  at  Purdue  University  on  fabric-reinforced  sands  (e.g., 
Reference  5),  and  its  properties  were  well  known.  The  grain  size  distribution 
of  the  sand  is  shown  in  rigure  19.  Other  classification  parameters  are: 


.gure  17.  Loading  Frame  (Reference  47). 
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General  View  of  Test  Box,  Data  Acquisition  System, 
and  Controls  for  t.he  MTS  Loading  System. 
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U.S.  Standard  Slava  Numbers 


Figure  19.  Grain  Size  Distribution  of  Ottawa  Flint  Sand 


1.6 
u.94 
0 . 34  nm 
2.65 

114 -9  Pcf;  er,n‘  n  =  J-44 
93.3  pcf;  o  r  U./7 

K  hid* 

The  grain  size  distribution  shown  in  t'i^re  19  is  tne  average  ot  .2 
separate  determinations,  which  varied  only  about  1  percent  at  each  sieve  mcc. 
these  oetermi nat 1 ons  were  made  in  order  to  verity  that  the  sand  stored  in 
several  containers  in  the  laboratory  was  the  same  as  that  used  previously.  '-Iso 
snown  in  Figure  19  is  the  gradation  supplied  oy  the  Ottawa  Silica  to. 

The  angle  of  internal  f notion  or  the  sand  was  determined  t run  eotr  ■  * 

shear  and  tri axial  compression  tests,  direct  shear  test  results  are  sr  ;w  - 

;  i  gure  23.  The  tests  were  multistage;  that  is.  only  one  test  spc«-iner  **a 
up  tor  each  run.  When  the  peak  of  the  shear  stress  oef  on  rat  i  or-  curve  *  t 
reached,  the  next  increment  of  normal  stress  was  applied.  A-  shown  ; r  •  • 
the  failure  envelopes  are  curved,  as  expected  for  sands  under  1  •* 
stresses.  Multistage  tests  may  also  contribute  to  envelope  curvy  „r<  , 
direct  shear  test  specimens  were  initially  set  up  at  a  dry  dens  it  v  ,  r  .  . 

which  is  very  close  to  the  average  density  cf  the  sand  u-  th*.  test  .  • 
pcf  or  a  relative  density  of  about  7u  percent) . 

The  normal  stress  existing  at  trie  elevation  ot  tut  tape  • 
specimens  is  only  on  the  order  of  iu  to  :•>'  psf ,  or  si  gm  t  i  r  i"t  1  v  V  •  . 
lowest  normal  pressure  at  which  tno  direct  shear  tests  «e*v  ts.t'  ' 

Attempts  to  Obtain  friction  angles  at  oven  lower  normal  -ft  •  •  «i  r,  . 

r.«l.  apparently  due  to  the  "internal"  friction  ot  V<:  shear  m.i.,  r--  ne  if  * 
precision  direct  shear  apparatus  (f.arol -warner)  was  used  for  t.  •  .. 

the  sliding  parts  were  specially  remachined  to  make  the  as  smooth  at  ....  • 

TnLn  a  teflon  -pased  antifriction  material  was  sprayed  on  c  c.n  mow  .  , 

t  •  ;'*mal  stress  of  leu  psf,  the  friction  of  the  apparatus  itself  was 
o  ■■>* .  but  tne  net  ooserved  shear  resistance  of  the  sand  machine  ‘n.to;  < 
traded)  ranged  from  2°  psf  to  145  psf.  The  two  data  points  shewn  on  r  i  < 

are  average  values.  Additional  direct  shear  tests  at  normal  stresses  '  .-.t  ■  *• 

j.  psf  were  not  attempted  because  of  the  undesirably  large  experiment  a1 
d  I*  ter. 

If  the  extrapolation  of  the  tfohr  failure  envelopes  to  the  origin  >  •  a? 

:  S’t.d  lines  in  Figure  2o )  is  correct,  then  the  friction  angle  at  the  re  ‘  ■  ve 
d-d  it/  of  7U  percent  for  the  confining  pressures  operative  in  the  test  s. 
wca'.i  approach  45  degrees.  At  higher  normal  pressures,  the  average  frieti 
,i"..l e  was  around  35  degrees,  a  value  more  typical  for  the  poorly  graded  rounoeci 
fra' n  Ottawa  sand.  The  direct  shear  results  in  figure  2d*  at  normal  stresses 
yo  iter  than  2u0  psf*  were  verified  by  tri axial  teste  conducted  earlier  on  this 
>j'c‘  s0nd  as  part  cf  ongoing  academic  work  at  i’urdue  .Ini  versit/ .  From  pi  jfe 
texts  on  unrei n forced  sands  described  later  in  this  report.,  it  is  possible  to 
f  :  1 ;  .ilate  1  friction  angle  from  clu'Sical  bearing  capacity  theory.  2s  i  ng  the 
r-  *  caghi -■■Ywer'hi  ft  bearing  capacity  factors  and  appropriate  empirical  correc¬ 
ts  ns  for  tne  sr.ape  of  the  footing,  an  average  fn,  fen  ingle  ,-f  42  degrees  w  is 
■  il.ijf  ated  tor  all  three  loaning  plates  tested.  These  r<  .u  1 1  >  con  f  i  rm  that,  tne 
;  •  : ;  *  i  s  w  i .  •  )  ri-i»  •  -«i  .n*  •  i«*  i  y  ♦  it  *■  r.  *  t  •  ”  ir  -  *  ‘  i  . 
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d.  Reinforcement 


Two  types  of  rei  nforcement  were  used  in  this  study.  The-  first  was  a 
hi gh-strength ,  high-modulus,  woven  polyester  fabric  which  has  been  usea  in  the 
past  for  other  studies  at  Purdue  University  or  yeotextiies  (see  References  A,  H, 
/,  and  4b).  The  fabric  is  typical  of  woven  mu  I  ti  f  i  1  ament  geotextilta  ;  :  o- 

duced,  for  example,  by  hi  col  on,  Carthage  Mills,  etc.,  although  this  fabric  actu¬ 
ally  was  obtained  from  AH  fodervavnarier ,  of  boras,  Sweden).  The  fabric  is  an 
industrial  grade,  woven  multi  filament  iiianuf  actured  from  polyester  fibers 
obtained  from  either  Hoechst  (Trevira'  No.  710 )  or  I Cl  ( " Terel in' " ) .  The  basic 
fiber  is  a  Hub  dtex  (Q90-ilenier)  polyester  with  a  20/ZbU  yarn  structure,  hide 
strip  tensile  tests  conducted  as  pare  of  this  study  gave  an  ultimate  tensile 
strength  of  about  lbC  Its/inch  Typical  tensile  lead-strain  curves  are  shown  in 
Figure  ?l.  In  this  fig.-re,  strain  was  determined  in  two  ways,  ( i )  as  the  ratio 
cf  the  crosshead  movement  to  the  crosshead  spacing  before  the  start  of  loading; 
or  12)  from  frequent  measurements  with  a  micrometer  of  two  marks  originally  2.5 
inches  apart  near  the-  center  of  the  test  specimen.  If  crusshead  strain  is  used 
as  the  criterion,  then  elongation  at  tot 5 ore  is  about  12  percent.  The  secant 
tensile  modulus  at  10- percent  strain  is  about  800  pounds/ inch.  It  should  be 
noted  that  the  data  shown  in  figure  21  were  obtained  from  a  modified  wide  strip 
tensile  test.  Although  the  test  is  not  yet  an  official  ASTM  standard,  it  is 
Currently  under  >wiew  by  AglM  '  muni t tee  013. cl  bib. 19  on  Geotextiles.  Approval 
is  expected  shortly,  'he  draft,  standard  calls  for  a  specimen  width  of  <s  inches, 
"est  specimens  were  ?  inches  wide.  This  difference  is  considered  negligible  for 
a  w-;ven  fabric.  From  previous  studies  (Rei ore nee  4),  the  stress-strain  behavior 
;s  almost  identical  for  both  the  warp  and  filling  directions. 

In  addition  to  high  tensile  strength  and  favorable  modulus,  the  fabric  nas 
several  other  important  properties.  ' t  has  excellent  creep  resistance.  Tests 
conducted  in  Sweden  (Retereno-  4s  j  .  . ndicatc-a  that  the  extrapolated  creep 
strength  after  more  than  l  .year  was  about  bu  percent  of  its  short-term 
strength.  After  1  year,  the  creep  was  fr.und  to  increase  about  0.18  percent  per 
log  cycle  of  time.  The  fiber  has  excellent  resistance  to  aging,  sunlight, 
weathering,  rotting,  bacteria,  and  rodents.  It  is  very  resistant  to  acids, 
relatively  resistant  to  bases  and  is,  in  general,  insoluble  in  most  inorganic 
solvents,  including  jet  and  diesel  fuels. 

The  other  reinforcirg  materials  used  in  this  study  were  “yeognds".  asi- 
(,-jlly  geogrids  are  extruaed  plastic  sheets  which  have  high  tensile  strength  and 
modulus.  They  are  made  from  a  high-density  polyethylene  or  polypropylene  in 
which  the  polymer  is  strongly  oriented  to  obtain  high  tensile  strengths.  On  a 
weight  basis,  the  grids  are  stronger  than  steel.  The  bas’c  polymers  are  highly 
resistant  to  chemical,  biological,  and  ultraviolet  radiation.  The  geogrids  used 
in  this  research  are  manufactured  in  England  by  the  Notion  Corporation  under  the 
trade  name  'ensar  .  Samples  were  obtained  from  Gulf  Canada,  Ltd.,  who  is  the 
current  North  American  licensee  and  who  will  scon  t>e  manuf ucturi ng  the  grids  in 
Canada,  data  provided  ny  the  manufacturer  are  given  in  Table  3  and  Figures  22, 
La  md  24.  For  reference  purposes,  physical  and  mechanical  data  for  all  or  the 
aval ’able  grids  are  shown,  although  only  types  SR?  and  552  were  tested.  The  ten¬ 
sile  strength  of  SR?  is  slightly  more  than  450  pounds- inch  at  a  maximum  straii 
of  12  percent  compared  to  lob  pounds/inch  for  the  woven  polyester.  The  tensile 
strength  of  SS2  is  somewhat  greater  than  that  of  the  polyester,  or  20b 
pounds' i nch  at  about  toe  same  ultimate  strain.  For  the  552,  note  that  the 


Load,  Ibf/in. 


TABLE  3.  PHYSICAL.  AND  CHEMICAL  PROPERTIES  OF  TENSAR*' 
GEOGRIOS 


A.  Polymer  Characteristics 

SRI  and  SR?  SSi  ana  SS2 

Polymer  High-density  Polypropylene 

polyethylene 

Snore  hardness  Lj 

(Din  53505)  b7  74 

Vicat  softening  point 

(Din  53460) (C)  127  1 48 

Impact  Strength  . 

(Din  53453 ) ( kJ /m  )  13.2  4.5 

Abrasion  resistance 

(Din  53754E)  14.0 

{nrr/10U  revs)  10 

Chemical  resistance  Resistant  to  all  naturally 

occurring  alkaline  and  acid¬ 
ic  conditions 

Biological  resistance  Resistant  to  attack  by 

bacteria,  fungi  and  vermin 

Sunlight  resistance  Resistant  to  UV  attack 


TABI£  3. 


PHYSICAL  AND  CHEMICAL  PROPERTIES  OF  TENSAP 
CEOGPIDS  (CONTINUED) 


Mechanical  Properties 


SRI 

SR2 

Tensile  strength  - 
maximum  (kN/m) 

84.0 

79.0 

Extension  at  maximum 
lead  ( c ) 

12.3 

12.0 

Extension  at  40 C  load  (%) 

3.5 

3.0 

Modulus  in  tension 
( N/m  ) 

5 

9 

.2  x  10 

4. 

1  x  10^ 

Thermal  stability 

Stable  over 
of  -60° 

temperature  range 

to  80  °C 

SSI 

SS2 

Across 

Roll 

width 

Along 

Roll 

length 

Across 

Roll 

width 

Along 

Roll 

length 

"Characteristic" 
tensile  strength  per 
metre  width  (kN/m) 

20.9 

12.6 

32.0 

18.0 

(Samples,  3  junctions  long  and  1  rib  wide  were  extended  at  a 
constant  rate  of  5L'nm/min,  at  a  temperature  of  20-l°C. 


Figure  22.  Tensile  Stress-Strain  Relationship  for 
fa)  SRI  and  f b )  SR2  Geogrids. 


Figure  23.  Results  of  Creep  Tensile  Tests  on 
(a)  SRI  and  (b)  SR2  Geogrids. 


Shear  Street  kN/m’ 


strength  in  the  direction  of  the  roll  length  is  about  60  percent  of  the  polyes¬ 
ter,  or  about  97  pounds/inch. 

The  data  given  in  Figure  24  are  from  Reference  49,  and  they  compare  the 
frictional  resistance  of  Tensar®  grids  with  that  of  other  materials  including  a 
"woven  polyester  civil  engineering  fabric"  (Reference  4)  which  is  the  same 
fabric  used  in  the  present  research. 

Geogrids  could  be  called  “second  generation"  geotextiles,  because  they  are 
in  general  much  stronger  than  most  geotextiles,  yet  their  cost  is  about  the  same 
for  the  same  weight  per  unit  area.  Secause  of  the  large  openings  in  the  grids 
which  can  provide  considerable  int  flock  resistance,  geogrids  have  an  advantage 
in  terms  of  frictional  resistance  (Figure  24). 

e.  Loading  Plates 

The  loading  plates  used  were  rigid  steel  plates.  Rigid  plates  were 
chosen  because, at  aircraft  tire  pressures,  the  tire  appears  to  be  essentially 
rigid  to  the  soil.  The  two  circular  loading  plates  were  3  inches  and  b  inches 
in  diameter.  The  3-inch  plate  was  0.75  inches  thick  while  the  b-inch  plate  was 
1  inch  thick.  The  plane  strain  plate  was  3  inches  wide  by  30  inches  long  py 
1  inch  thick  and  was  stiffened  along  its  entire  length  to  minimize  the  possibil¬ 
ity  of  bending  during  the  test.  The  plates  were  attached  to  the  load  cell  with 
reducing  couplings  ana  tnreaded  studs. 

Photographs  of  the  loading  plates  will  be  shown  later. 

f.  Instrumentation  and  Data  Acquisition  System 

In  this  paragraph.only  a  brief  description  will  be  given  of  the  instru¬ 
mentation  utilized  in  the  test  program.  More  detailed  descriptions  will  be 
given  later  in  this  chapter  when  describing  the  test  procedure. 

The  deformation  of  the  load  plates  was  measured  either  by  an  external 
deformation  transducer  (DCDT,  direct  current  differential  transducer)  or  by  the 
LVDT  (linear  variable  differential  transformer)  in  the  MTS  system.  The  external 
instrument  was  a  Hewlett-Packard  Model  7DCDT  3.000,  with  a  linear  stroke  range 
of  -o.-jL  inches. 

The  lead  applied  by  the  actuator  ram  was  measured  either  by  a  Sensotec 
Model  41  load  cell  with  a  capacity  of  50,000  pounds  or  a  Lebow  Model  16.103  load 
cel!  with  a  5, 000-pound  capacity,  both  these  load  cells  are  strain  gage  type. 

Deflections  of  the  sand  surface  were  measured  by  several  DCDTs  mounted  in  a 
beam  bolder  and  running  radially  from  the  edge  of  the  loading  plate  to  one  side 
of  the  box.  All  OCDTs  were  powered  by  a  Hewlett-Packard  Model  62U5B  power  sup¬ 
ply.  Figure  25  shows  the  mounting  beam  for  the  DCDTs  as  well  as  the  load  cell 
and  plate  for  a  b-inch  CP  test.  The  3-inch  DCDT  can  be  seen  just  to  the  left  of 
the  load  cell.  As  shown  in  the  photograph,  the  DCDT  support  beam  was  held  by 
long  threaded  bolts  to  "box  girders”  supporting  it  to  facilitate  leveling  of  the 
DCDi  support  beam. 


izure  25.  DCDT  Support  Beam,  Load  Cell,  and  6-inch  Diameter 
Load  Plate. 
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The  location  (radial  distance)  of  the  DCDTs  for  each  test  varied  somewhat, 
depending  on  the  test  series  and  the  size  and  geometry  of  the  loading  plate. 
Figure  26  and  Table  4  show  the  plan  and  location  of  the  DCDT  for  each  test 
series. 

An  Analog  Devices  Inc.  MACSYM  2  data  acquisition  system  (DAS)  was  used  to 
facilitate  rapid  data  acquisition  and  reduction.  The  system  is  a  "smart  data 
logger."  It  uses  software  (a  variation  of  the  computer  language  BASIC)  to  both 
acquire  and  operate  on  the  data.  To  obtain  hard  copy  output  from  the  MACSYM,  a 
Heathkit  printer,  Model  H- 14  was  connected  to  the  DAS.  Figure  27  is  a  photo¬ 
graph  of  the  data  acquisition  system  and  printer  unit.  Figure  28  shows  a 
schematic  diagram  for  connecting  the  instrumentation  to  the  MACSYM. 

Additional  instrumentation  included  a  number  of  Micromeasurement  type  large 
strain  SR4  strain  gages.  It  was  originally  intended  that  the  strain  gages  would 
be  attached  to  the  fabric  at  several  locations  on  the  test  specimen.  However, 
considerable  technical  difficulties  arose  during  some  initial  pilot  tests  with 
the  strain  gages  attached  to  fabric  tensile  test  specimens.  Because  tne  test 
results  appeared  so  erratic  and  unreliable,  it  was  decided  not  to  pursue  this 
line  of  research  further.  Tests  results  indicated  that  the  stress-strain  pro¬ 
perties  of  the  fabric  were  altered  significantly  by  the  presence  of  the  gage  on 
the  fabric,  probably  because  of  the  epoxy  type  adhesive  used.  (Coincidentally, 
the  problem  of  large  strain  measurements  on  geotextiles  is  considered  a  No.  1 
Priority  Research  Need  by  the  Committee  on  Soil  and  Rock  Instrumentation  of  the 
Transportation  Research  Board.)  Unfortunately,  we  were  not  able  to  make  a  sig¬ 
nificant  step  toward  solving  that  problem. 

A  Bison  strain  indicator  and  soil  strain  gage  proved  unsatisfactory  within 
the  constraints  of  the  study.  Some  difficulty  was  encountered  in  calibrating 
these  strain  gages  until  very  late  in  the  research.  Furthermore,  it  was  not 
technically  possible  to  read  several  soil  strain  gages  at  once  during  a  test  and 
avoid  the  interference  between  gages  in  the  test  box. 

The  only  other  measurement  made  during  the  test  was  the  movement  of  the 
side  wall  of  the  box.  This  deformation  was  measured  by  an  ordinary  dial  indica¬ 
tor,  accurate  to  0.01  millimeters.  These  measurements  will  be  reported  later  in 
this  chapter. 


4.  FABRICATION  AND  CONSTRUCTION  OF  TEST  SPECIMENS 
a.  Sand  Placement  and  Density  Control 

There  was  considerable  initial  concern  about  being  able  to  place  the 
large  volume  of  sand  in  the  test  box  at  a  reasonably  uniform  density.  It  was 
desirable  that  the  sand  be  on  the  dense  side,  that  is,  not  loose.  Several  sys¬ 
tems  were  considered  for  both  sand  placement  and  density  control .  It  was  ori¬ 
ginally  anticipated  that  the  sand  could  be  rained  in  the  box,  and  by  keeping  the 
height  of  fall  constant  as  the  thickness  of  the  sand  built  up,  the  resulting 
density  could  be  maintained  constant.  Also  considered  was  the  air-activated 
sand  spreader  system  described  by  Butterfield  and  Andrawes  (Reference  50),  and 
the  scheme  used  previously  at  Purdue  University  by  Brummond  and  Leonards  (Refer¬ 
ence  51).  All  these  schemes  were  rejected  because  of  the  presence  of  the 
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Figure  26.  Location  of  DCDTs  for  Each  Test 
Configuration  (see  also  Table  4). 


TABLE  4.  LOCATION  OF  DCDT  FOR  EACH  TEST  CONFIGURATION 
(SEE  ALSO  FIGURE  2b) 
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Figure  28.  Schematic  Diagram  of  Test  Instrumentation 
and  Data  Acquisition  System. 


reaction  frame  (Figure  16).  Either  the  reaction  frame  would  have  to  be  moved 
before  each  test,  or  the  loading  box,  once  it  was  filled  with  sand,  would  be 
moved  back  onto  the  base  of  the  loading  frame  for  testing.  Both  schemes  were 
deemed  impractical.  The  weight  of  the  reaction  frame  was  several  hundred  pounds 
and  was  securely  bolted  to  the  base  plate.  It  seemed  impractical  to  repeatedly 
bolt  and  unbolt  the  frame  prior  to  a  test.  The  weight  of  the  sand  plus  box  was 
in  excess  of  1200  pounds.  Consequently,  the  placement  scheme  which  was  finally 
adopted  is  as  follows: 

The  sand  was  placed  in  three  layers,  each  approximately  6  inches  thick  and 
vibrated  after  each  layer.  The  vibration  was  carried  out  with  an  FMC  Syntron 
Magnetic  Vibrator  Model  V51C1,  controlled  by  a  Syntron  Electronic  Controller 
Model  SCR-18 .  The  vibrator  was  attached  to  a  20- i nch  by  20-inch  by  1/4-inch  thick 
steel  plate.  Figure  29  shows  the  plate  vibrator  in  position  prior  to  compac¬ 
tion  of  the  first  layer  of  sand.  The  ropes  which  were  attached  to  the  rubber 
pads  were  used  to  facilitate  moving  the  vibrator  and  plate  to  the  four  positions 
in  the  box.  The  frequency  of  vibration  was  60  Hert2.  The  control  was  set  on 
"maximum"  at  all  times,  and  the  measured  amplitude  of  the  plate  vibration  was 
0.056  inch.  The  vibration  pattern  is  shown  in  Figure  30.  The  area  of  the 
vibratory  plate  was  such  that  there  was  about  a  9-inch  overlap  in  the  center 
section.  The  initial  loose  lift  heights  were  approximately  6.5  inches.  After 
vibration,  each  layer  ended  up  about  6  inches  thick.  Thickness  control  was 
maintained  by  a  scale  on  the  inside  of  the  box. 

The  same  pattern  of  vibration  was  followed  for  the  tests  in  which  rein¬ 

forcement  was  used.  In  these  cases,  however,  an  extra  pass  of  vibration  was 
applied  to  the  layer  of  sand  on  top  of  the  fabric.  As  will  be  shown  later, 
average  densities  of  the  test  where  fabric  reinforcement  was  used  tended  to  be 
slightly  greater  than  the  densities  of  the  tests  without  fabric  rei nforcement, 
probably  due  to  this  extra  time  of  vibration. 

Oensity  measurement  and  control  were  originally  attempted  with  a  nuclear 
density  gage,  but  the  readings  were  too  erratic,  probably  due  to  all  the  high 
density  material  (steel)  in  the  area  of  the  test  box  and  load  reaction  frame. 

To  obtain  accurate  results  with  nuclear  density  meters,  no  other  materials  other 
than  the  specimen  under  test  should  be  in  the  vicinity. 

The  two  approaches  that  were  finally  used  were  to  measure  (1)  the  average 

overall  density  of  the  sand  in  the  box,  and  (2)  the  density  of  individual  small 
cans  placed  in  several  locations  in  the  box.  First,  the  overall  density  was 
obtained  by  carefully  weighing  all  the  sand  placed  in  the  box.  This  procedure 
was  done  twice  and  the  average  density  was  found  to  be  107.8  pounds/cubic  foot. 
This  is  equivalent  to  a  relative  density  of  slightly  more  than  70  percent.  The 
second  procedure  was  followed  in  23  tests.  Small  metal  (water  content)  cans 
were  placed  in  several  locations  in  the  test  box  (Figure  31).  To  determine  the 
volume  of  these  density  cans,  each  was  filled  with  distilled  water  of  known  tem¬ 
perature  and  the  can  was  carefully  weighed.  The  measurement  was  repeated  five 
or  six  times  and  the  average  volume  used.  All  weights  were  determined  on  a  pre¬ 
cision  laboratory  balance.  Placement  of  the  cans  in  the  box  was  according  to 
the  following  scheme: 


Figure  30.  Vibration  Pattern. 
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*  Exception:  for  PS  tests;  this  can  was 
positioned  near  the  south  wall. 

The  cans  were  placed  aoout  4  inches  from  the  wall*  except  for  the  can  under  the 
centerl i ne . 

At  the  conclusion  of  the  loading  test,  the  cans  were  carefully  excavated, 
struck  off  level,  weighed,  and  the  density  of  the  sand  in  each  can  was  calcu¬ 
lated.  The  results  of  these  measurements  are  given  in  Table  5.  The  average 
value  of  all  the  measurements  by  this  scheme  was  slightly  less  than  the  overall 
average  density  in  the  box  of  107.9  pcf.  The  range  of  densities  measured  in 
each  test  was  between  about  1  and  4  pcf.  There  seemed  to  be  no  definitive  pat¬ 
tern  of  the  densi ties, al though ,  with  a  few  exceptions,  the  two  bottom  cans 
tended  to  be  lower  than  those  on  the  middle  and  the  top  layers.  Occasionally, 
the  middle  layer  cans  appeared  to  have  the  higher  density  (as  might  be  expected 
from  work  with  field  vibratory  rollers).  The  surface  of  the  sand  after  compac¬ 
tion  appeared  to  be  very  dense  to  the  touch. 

After  the  last  layer  was  compacted  on  tne  rei nforcement,  a  smal 1- di ameter 
wire  probe  was  used  to  verify  the  depth  to  the  fabric  layer. 

At  the  conclusion  of  a  test,  all  but  2  to  4  inches  of  the  sand  was  removed 
and  placed  in  storage  cans.  The  sand  left  in  the  bottom  of  the  box  was  stirred 
repeatedly.  For  preparation  of  the  next  test,  the  sand  from  storage  was  loosely 
dumped  into  the  box  so  that  the  initial  loose-lift  thickness  was  around  b.5 
inches.  Then  the  compaction  process  proceeded  as  described  above  in  preparation 
for  the  next  test. 

b.  Fabric  and  Geogrid  Placement 

The  fabric  or  geogrids  were  placed  at  depths  and  in  layers  as  described 
previously  (see  Table  1  for  the  test  code).  The  rei nforcement  was  placed  on  the 
compacted  sand  surface  and  smoothed  by  hand  so  that  it  was  in  full  contact  with 
the  densified  surface.  Figure  32  shows  the  sand  being  loosely  dumped  on  top  of 
a  fabric  layer.  Then  this  final  layer  of  sand  was  compacted  as  described  above. 
In  the  majority  of  the  tests,  the  fabric  or  geogrid  reinforcement  was  not 
attached  to  the  walls  or  pretensioned  as  was  done  by  Hal i burton  ,  Lawmaster  and 
King  (Reference  13).  It  was  fe  t  that  tests  of  the  loosely  placed  fabric  under 
the  low  confining  pressures  applied  by  the  thin  layer  of  sand  would  provide  a 
lower  bound  for  the  effect  of  the  rei nforcement.  For  comparison  purposes,  how¬ 
ever,  a” few  tests  were  carried  out  with  (1)  partial  fixity,  and  (2)  total  fix¬ 
ity.  The  results  of  these  tests  will  be  described  later. 


TABLE  5.  DRV  DENSITIES  AS  DETERMINED  BY  THE  SMALL 
DENSITY  CAN  MEASUREMENTS 
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c.  Mounting  of  the  Load  Plate  and  Load  Cell 

After  the  fabric  and  final  layer  were  placed  and  densified,  the  load 
ceM  and  appropriate  loading  plate  were  attached  to  the  actuator  piston  by 
threaded  studs.  The  cable  from  the  load  cell  was  then  attached  to  the  MTS  sys¬ 
tem.  Figure  33  is  a  closeup  of  the  load  cell  and  6-inch  diameter  load  plate. 

d.  Instrumentation 

As  described  before,  the  DCOTs  were  set  up  according  to  the  pattern 
shown  previously  in  Figure  26  and  Table  4.  In  the  closeup  view  of  Fiyure  33, 
the  3-inch  DCDT  is  positioned  just  prior  to  loading  of  the  6CP  tests  (Series  II 
of  Figure  26[b]).  Figure  34  is  a  closeup  of  the  DCDT  support  beam.  The  cores 
of  the  DCDT  rest  directly  on  small  cardboard  discs  on  the  surface  of  the  sand  to 
record  the  upward  and  downward  movements  of  the  sand  surface  as  the  loading 
plate  is  pressed  into  the  surface. 

After  all  the  instrumentation  was  mounted  and  appropriate  connections  made 
to  the  MTS  and  the  data  acquisition  system,  the  positions  of  the  DCDTs  were 
adjusted  so  that  the  cores  were  within  the  linear  range  of  the  expected  movement 
during  the  test.  A  MAC3YM  program  was  written  to  assist  in  this  positioning. 

All  connections  were  checked  to  make  sure  everything  was  ready  for  the 
loading  test.  The  DAS  was  turned  on,  and,  after  a  brief  warmup  period,  the 
piston  of  the  MTS  was  lowered  to  apply  a  very  small  seating  lead  (1  to  2  pounds; 
on  the  surface  of  the  sand.  Now  the  test  was  ready  to  begin. 

e.  Loading  Test 

Loading  was  accomplished  with  the  MTS  system  in  "load  mode."  The 
hydraulic  pressure  range  was  set  so  that  sufficient  load  can  be  applied  by  the 
actuator.  The  load  range  on  the  MTS  was  set  to  an  estimated  value  of  the  max¬ 
imum  load  so  that  maximum  sensitivity  of  the  load  cell  can  be  achieved.  The 
load  cell  and  the  MTS  LDVT  were  "zeroed"  as  close  as  possible  with  digital 
voltmeters  (DVM,  Figure  28).  To  apply  the  load,  the  Set  point"  on  the  MTS  Ser- 
vac  control  unit  was  turned  slowly.  The  DAS  was  programmed  to  take  readings  of 
all  channels  at  a  specified  frequency.  For  the  present  tests,  readings  were 
made  every  0.25  or  0.5  seconds.  Thus,  for  each  test,  about  1000  readings  were 
taken  of  all  channels  and  stored  in  the  memory  of  the  DAS.  The  load  was  applied 
at  approximately  a  constant  rate  by  monitoring  the  output  from  the  load  cell  on 
the  DVM.  The  load  was  applied  so  that  failure  was  reached  in  7  to  8  minutes. 

Failure  occurred  suddenly  in  most  cases,  and  sometimes  additional  loao  was 
applied  to  see  if  there  was  a  work-hardening  effect.  At  the  end  of  each  test, 
the  loading  piston  was  retracted,  the  plate  raised,  instrumentation  removed,  and 
photos  and  or  sketches  made  as  appropriate.  After  excavation,  appearance  of  the 
fabric  and/or  grids  was  noted.  For  the  fabric  tests,  there  was  no  obvious  phy¬ 
sical  evidence  of  disturbance  or  abrasion  on  any  test  specimens  on  the  fabric. 

On  the  plane  strain  specimens  on  the  yeoyrids  (Tensar ”  SR2 ) ,  a  definite  "per¬ 
manent  set"  of  a  "dip"  directly  under  the  loading  plate  was  shown.  This  can  De 
seen  in  Figures  35  and  36. 


During  the  test,  observations  were  made  of  the  deformation  of  the  side  of 
the  box  by  means  of  a  mechanical  dial  indicator  accurate  to  the  nearest  0.01 
millimeter.  These  measurements  will  be  discussed  in  the  next  section. 

5.  EXPERIMENTAL  RESULTS 

The  primary  test  variables  were  described  earlier  in  this  chapter  (see 
Tables  1  and  2).  It  was  originally  planned  to  conduct  tests  with  vibratory  and 
impact  as  well  as  quasi-static  type  loadings.  The  analytical  work  (Section  II) 
did  not  require  other  than  quasi-static  loading;  consequently,  other  modes  of 
loading  were  not  utilized  in  these  series  of  tests. 

a.  Photographs  of  Tests  After  Failure 

Figure  37  shows  a  test  with  fabric  rei nforcement  using  a  3-inch  diame¬ 
ter  loaded  plate  after  failure.  Figure  38  is  a  closeup  view  of  the  same  test. 
Note  the  uniform  and  symmetrical  appearance  of  the  failure  surface.  The  next 
photograph.  Figure  39,  shows  the  failure  surface  after  the  piston  has  been  with¬ 
drawn.  Again,  note  the  uniformity  of  the  surface.  These  observations  tend  to 
verify  that  the  sand  density  in  the  vicinity  of  the  piston  was  sensibly  uniform. 

The  next  three  photographs  show  plane  strain  tests.  Figure  40  shows  the 
plane  strain  test  setup  prior  to  loading.  The  next  photograph  (Figure  41)  shows 
the  same  test  after  failure.  At  the  "north"  side  of  the  box,  that  is,  under  the 
far  end  of  the  DCDT  support  beam,  it  can  be  seen  that  the  fabric  has  pulled  away 
from  the  edge  of  the  box,  leaving  a  depression  in  the  sand  at  the  edge.  This 
was  typical  behavior  of  both  the  circular  and  plane  strain  tests  at  failure. 

The  amount  of  "pull"  away  from  the  edge  of  the  box  depended  on  the  total  defor¬ 
mation  of  the  loading  plate.  Figure  42  shows  a  plane  strain  test  with  a  view 
from  above.  Note  that  the  beam  appears  to  be  slightly  off-center,  an  occurrence 
typical  of  the  plane  strain  tests.  The  beam  was  set  perpendicular  to  the  sides 
of  the  box  at  the  start  of  the  test.  When  failure  occurred,  some  slight  rota¬ 
tion  of  the  beam,  either  to  the  left  or  to  the  right,  was  usually  observed. 

For  the  partial  and  full  fixity  tests,  the  geotextile  was  held  by  means  of 
wood  slats  bolted  to  the  box  sides.  Figure  43  shows  the  fabric  clamped  in  place 
prior  to  placement  of  the  final  sand  layer.  Figure  44  shows  how  the  sand  and 
fabric  have  pulled  away  from  the  right  side  wall  and  clamp  in  a  circular  plate 
test  with  partial  fixity.  In  Figure  45,  after  failure,  the  sand  has  been 
removed,  and  an  arrow  drawn  on  the  fabric  to  indicate  how  much  the  fabric  has 
moved  (in  this  case,  about  1  inch).  Again,  this  was  a  test  with  the  fabric  par¬ 
tially  fixed;  that  is,  the  fabric  was  folded  under  about  1  inch  all  around  prior 
to  placement  of  the  final  sand  layer.  No  photographs  were  taken  of  the  test 
with  full  edge  fixity  because  its  appearance  was  so  similar  to  the  tests  with 
partial  fixity. 

b.  Load  vs.  Settlement  of  Plates 

Test  results  are  shown  in  the  form  of  the  load-settlement  response  of 
the  loading  plate  size  and  shape,  both  with  and  without  rei nforcement.  The 
effect  of  the  several  variables  on  the  test  resul  will  now  be  discussed. 
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Figure  46  shows  the  load  vs.  settlement  behavior  for  the  3- inch  diameter 
and  the  6- inch  diameter  plate  tests. 

The  increase  in  the  bearing  capacity  due  to  the  larger  plate  size  is  obvi¬ 
ous.  Note,  too,  the  increase  in  bearing  capacity  due  to  the  presence  of  the 
reinforcement.  The  amount  of  this  increase  as  a  function  of  depth,  type  and 
number  of  layers  will  be  discussed  later.  Also,  there  was  a  significant  work¬ 
hardening  effect  in  the  reinforced  sands  as  the  load  continued  to  be  applied  to 
the  plate  after  first  yielding  had  occurred.  In  almost  every  test,  a  fall-off 
of  load  was  observed  as  deformation  continued;  then,  if  the  test  was  carried  far 
enough,  the  reinforced  soil  system  began  to  take  up  load  again.  In  most  cases, 
the  range  of  the  actuationLVDT  was  exceeded,  so  that  measured  settlements  were 
not  obtained  beyond  about  1.5  to  2.0  inches.  The  work-hardening  effect  was  not 
so  apparent  for  the  unreinforced  sands,  although  a  peak  in  the  load-settlement 
curves  can  be  seen.  Such  behavior  is  not  unexpected  for  dense  sands. 

The  load-settlement  results  of  Hal i burton,  et  al .  (References  12  and  13) 
indicated  a  much  greater  work-hardening  effect  than  the  present  tests  show, 
probably  because  of  differences  in  the  edge  fixity  conditions.  This  point  will 
be  further  discussed  later. 

It  is  encouraging  from  Figure  46  to  note  that  the  test  results  are  reason¬ 
ably  repeatable,  both  with  respect  to  modulus  and  to  maximum  load.  This  effect 
is  apparent  in  Figure  47,  where  the  deflection  scale  (ordinate)  is  very  much 
expanded.  The  initial  moduli  (or  more  properly  the  coefficient  of  subgrade 
reaction,  which  miyht  be  determined  from  a  plate  loading  test  similar  to  those 
carried  out  here)  as  determined  from  the  slopes  of  the  load-settlement  curves, 
are  approximately  the  same.  In  Figure  48,  where  the  effect  of  different  types 
of  reinforcement  is  shown,  the  initial  moduli  of  the  tests  with  the  woven  fabric 
reinforcement  are  seen  to  be  about  the  same  as  their  counterpart  tests  with  the 
geogrids.  The  moduli  for  the  2- inch  depth  of  burial  are  somewhat  less  than  those 
for  1-inch  depth  of  burial,  independently, of  the  type  of  the  reinforcement.  The 
presence  of  the  reinforcement  at  1-inch  depth  clearly  affected  the  initial 
modulus  (or  coefficient  of  subgrade  reaction). 

Figure  49  shows  the  effect  of  the  number  of  layers  of  the  fabric  reinforce¬ 
ment.  Tests  on  multiple  layers  were  only  conducted  using  the  6-inch  CP.  The 
difference  in  behavior  between  the  two  tests  with  two  layers  of  reinforcement  is 
quite  apparent,  but  it  is  not  obvious  why  that  difference  occurred.  Because  of 
time  and  other  constraints,  these  experiments  were  not  repeated;  therefore,  it 
is  not  known  whether  the  effect  is  a  test  artifact  or  not.  The  work-hardening 
effect  is  not  so  apparent,  probably  because  of  lack  of  edge  fixity,  although 
some  effect  was  shown  on  the  two-layer  test  with  fabric  at  1.6  inerts  spacing. 

As  expected,  there  was  a  significant  increase  in  the  load-carrying  capacity  of 
the  two-layer  system  at  a  shallow  depth  of  burial  (the  6CP2W-G.8  test). 

Figure  50  shows  the  load-settlement  results  for  the  plane  strain  tests. 

There  was  a  work-hardening  effect  for  both  tests  with  the  fabric  reinforcement 
as  well  as  with  the  geogrid  reinforcement  at  1-inch  depth.  There  was  not  much 
difference  in  response  between  the  two  kinds  of  rei nforcement,  although  both 
showed  a  significant  improvement  over  unreinforced  sand.  The  moduli  of  all 
tests  with  reinforcing  were  slightly  greater  than  the  unreinforced  sand,  although 
the  difference  is  small. 
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Figure  51  shows  the  effect  of  edge  fixity.  The  second  run  of  test  3CPU- 
1- inch  I’F  had  a  disturbed  surface  and  consequently  failed  at  a  significantly  ler 
ultimate  load  than  the  first  run  or  that  test  i 1 02  pounds  vs.  56  pounds). 
Surprisingly,  there  was  very  little  di fference  between  the  tests  with  full  fix¬ 
ity  and  those  with  partial  fixity.  Recall  that  in  the  full  • i x i ty  test,  the 
geotextile  was  attached  by  wooden  slats  bolted  to  the  side  of  the  box.  For  tne 
partia1  fixity  test,  a  fold  of  about  1  inch  was  made  in  the  fabric  all  around 
tne  periphery.  The  full  fixity  test  was  expected  to  be  an  upper  bound  of  tnc 
effect  of  reinforcement,  while  the  test  with  partial  fixity  would  show  the 
i ntermediate  value  with  respect  to  full  fixity  and  no  edge  fixity.  These 
results  will  be  discussed  in  greater  detail  later  in  this  section. 

c.  Hiniate  Load  and  Stress 

Table  6,  an  expansion  of  Table  2,  describes  the  various  tests,  as  well 
as  some  test  results. 


'able  o  gives  ultimate  load  P  at  initial  yield,  P  .  and  at  a  plate  set¬ 
tlement  of  1  inch  for  each  test.  The  depths  to  the  reinforcement  layers  are 
shew"  as  a  ratio  of  the  plat“  diameter  (or  plate  width  in  the  case  of  the  plate 
strain  test),  or  a,b.  Tne  unit  bearing  stress  at  failure,  q,  is  P  /A.  From 
the  values  in  Table  t,  it  is  apparent  that  reinforcement  results  in  a  signifi¬ 
cant  increase  in  q  for  a  1 1  plate  sizes.  Significant  improvement  also  results 
with  the  addition  of  a  second  fabric  layer,  as  shown  in  the  bCPzW  tests.  fjll 
or  partial  en.,0  fixity  apparently  nas  vc-y  little  effect,  and,  surprisingly,  the 
use  of  tne  richer  strength  qeoyrids  did  not  increase  tne  ultimate  stress'd*  much 
as  x  xpecteu. 

:  inquet  ana  Lot-  (Reference  22)  expressed  the  degree  of  improvement  due  to 
reinforcement  in  terms  of  a  peering  capacity  ratio  (BCR).  They  defined  the  BCR 
as  r,  with  reinforcement  divided  by  j,  without  reinforcement,  at  the  same  verti¬ 
cal  deformation.  1 1  the  present  test  series,  initial  yield  occurred  at  about 
the  same  settlement,  and  once  failure  occurred,  the  amount  of  vertical  deforma¬ 
ble  v.:ur*-ir,g  thereafter  does  not  appear  to  make  much  difference  in  calculated 
k R  values  1  see  the  load-settlement  curves.  Figures  46  to  51).  Therefore, 
present  results  were  put  in  terms  of  a  modi f i ed  BCR,  wherein  the  q / q  ratio  was 
deters. i  ted  f o>  the  maximum  observed  values  of  each.  The  modified  BCR  values  for 
the  van  pus  Vtr  with  respect  to  the  test  on  the  unreinforced  sand  tor  the 
va’'iev>  plate  sizes  and  shape i  are  also  snown  in  Table  b. 

It  m  uf  .nterpit  to  compare  tne  results  cf  the  tests  on  unreinforced  sec- 
tics  with  tne  t! cf  i  :a5  bearing  capacity.  This  was  done  using  Terzaghi- 
M eye'*h;.ff  bearing  capacity  factors  (References  52  and  63).  Table  7  shows  tne 
theoretical  ultimate  bearing  capacity  in  pounds  for  the  three  test  configura¬ 
tions  for  a  range  of  friction  angles.  Also  shown  in  the  table  are  the  observed, 
maxinu":  values  of  load  on  the  plate,  lor  all  of  the  tests,  empirical  correc¬ 
tions  tin-  t.ne  shape  cf  the  loaded  area  were  made. 
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The  .1  j reenie nt  of  the  theoretical  values  for  ,  ;  4k  with  oDserved  values  is 
excellent.  This  value  of  *  also  agrees  with  the  results  of  the  direct  shear 
tests  presentee  in  Figure  CO  arid  verifies  the  observation  that  the  sands  were 
relatively  dense. 

ompa  r  i  son  of  present  results  for  reinforced  sands  with  previously  ;\.ib- 
'sned  theoretical  cevel opr.er.ts  is  not  easy,  dinquet  and  cee  (Referenc  -  .  k ) 
suggested  tnat  the  theoretical  bearing  capacity  ratio  for  one  reinforcement 
layer  should  be  about  Since  their  theory  was  for  a  strip  footing  and  for 

st-'ips  of  reinforcing  rather  tnan  sheets,  this  comparison  may  lack  validity, 
-’resent  L">  test  results  liable  o)  indicate  OCR  values  significantly  greater  than 
ninguet  and  Lee's  results  for  two  reinforcement  layers  suggest  a  theoreti¬ 
cal  LLP  of  about  l.S.  without  any  correction  for  strip  vs.  circular  footings, 
present  .•«.  A  values  tor  two  fabric  layers  were  significantly  greater.  1  he  only 
,tr-r  theoretical  model  for  hearing  capacity  of  fabric  reinforced  sands  has  been 
;  •>.  seated  Dy  rial  i  burton,  et  al  .  (Reference  13).  Their  results  will  be  discussed 
a  it  . 


■aliburton,  ct  al .  t Reference  Id  report  that  an  improvement  on  the  order 
in  fabric-soil  load  capacity  can  he  inferred  from  the  work  of  barenberg 
•■.ward,  and  tneir  coworkers  (References  kb,  kb,  t !  and  b4 ) .  Our  tests  verify 
cat  i  .>■  is  a  good  average  improvement  value. 

c.  effect  of  ’'enforcement  -c pin  and  spacirj 


F  i (jure  k  shows  the  maximum  lead  on  the  plates.  F  ,  vs.  the  depth,  n, 
*  rirst  re i n ‘or,. t ment  layer,  it  appears  that  depth  na^  very  little  iutlo- 
on  the  first  rei nf orcewent  layer  for  both  tne  ype  and  bC?’  tests  if  tne 
t  .  tr,  of  t-iiri al  is  less  than  inches,  for  the  btF  test,  when  the  depth  cachet! 

echos,  there  was  a  cijOiricant  decrease  in  the  k  •  On  the  ordinate  c  1  i  y- 
■v  Lk,  the  theoretical  bearing  capacity  values  form£fte  unreinforced  sand'  for  a 
c  4k  degrees  are  shown  r’Th.l 


T he  use  of  two  1  ay . r  .  of  woven  fabric  reinforcement  with  the  bit  tests 
ws  a  sign'ti  ant  iniprov»>r.*ort.  It  is  not  known  what  would  have  happened  if 
e**  rei  nf orcemcit  depths  had  boon  tested. 
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With  respect  to  increasing  the  depth  to  the  first  reinforcement  layer,  the 
results  for  the  circular  plates  are  exactly  the  opposite  as  with  the  plane 
strain  tests.  In  the  latter  tests,  the  increased  significantly  with  an 
increase  in  the  depth  of  the  reinforcement  layer,  at  least  to  the  depths  tested. 
It  is  likely  that  P  would  eventually  decrease  if  the  depth  to  the  first  layer 
was  sufficiently  laW)!.  The  exact  depth  at  which  this  decrease  would  begin  to 
occur  will  require  additional  testing. 

The  effect  of  the  geogrid  rei nforcement  for  both  the  CP  and  PS  tests  was 
surprising  in  view  of  the  differences  in  modulus  and  ultimate  tensile  strength 
Detween  the  grids  and  fabric.  Apparently  the  advantage  of  interlock,  which  the 
yrids  obviously  have,  was  not  realized  with  the  relatively  fine  Ottawa  sand 
l that  is,  "relatively  fine"  with  respect  to  the  dimension  of  the  "threads"  of 
the  grids).  With  the  fabric  reinforcement,  the  diameter  of  the  grains  is  of  the 
same  order  of  magnitude  as  the  diameter  of  the  threads,  and  therefore  signifi¬ 
cant  reinforcing  effects  occur  (References  4,  48,  49),  as  indicated  by  the 
marked  increase  in  bearing  capacity  of  the  reinforced  vs.  the  unreinforced  sands 
(Figure  52) . 

Similar  effects  are  also  shown  in  Figure  53,  where  the  BCR  is  plotted  vs. 
the  depth  in  inches  to  the  first  reinforcement  layer.  The  bearing  capacity 
ratio,  BCR,  is  defined  herein  as  the  bearing  capacity  of  the  reinforced  section 
to  the  bearing  capacity  of  the  unreinforced  section  (after  Reference  22).  Fig¬ 
ure  53  shows  that  the  BCR  is  much  greater  for  the  reinforcement  at  shallower 
depths  than  for  the  reinforcement  at  greater  depths.  In  both  the  3CP  and  the 
6CP  tests,  a  general  falling  off  of  the  BCR  occurs  as  the  depth  of  the  first 
layer  increases.  The  same  trend  is  apparent  for  tests  with  two  layers  of  rein¬ 
forcement  and  for  the  grid  reinforcement.  In  fact,  for  those  two  situations, 
the  decrease  of  BCR  with  increasing  depth  is  much  greater  than  for  the  single 
layer  of  fabric.  The  only  exception  to  these  trends  is  with  the  PS  tests,  and 
as  mentioned  above,  it  is  possible  there  would  be  a  decrease  in  BCR  if  the  depth 
to  the  first  layer  was  increased  sufficiently.  In  both  cases,  fabric  and  grids, 
the  BCR  increased  from  about  1.5  to  slightly  more  than  1.6  as  the  depth  was 
increased  from  1  to  2  inches.  When  the  effect  of  partial  and  full  fixity  of  the 
fabric  is  considered,  it  is  interesting  that  the  BCR  for  those  two  tests  was 
only  about  1.7  as  opposed  to  2.0  or  more  for  the  "free"  edge  test.  Why  this 
decrease  occurred  is  not  obvious,  but  from  the  analysis  of  Haliburton's  test 
data  (References  12  and  13)  presented  later  in  this  section,  it  appears  that 
edge  conditions  and  even  partial  fabric  prestress  do  not  significantly  influence 
the  bearing  capacity  at  initial  yield  of  the  loaded  plate. 

Also  shown  in  Figure  53  are  the  results  of  two  series  of  bCP  tests,  "old" 
and  "new"  series.  The  "old"  series  consisted  of  tests  which  were  conducted  with 
the  5U-kip  load  cell,  and,  consequently,  the  sensitivity  of  the  loads  measured 
left  something  to  be  desired.  Those  tests  with  fabric  at  1-,  2-  and  3-inch 
depths  were  repeated,  and  although  the  observed  strengths  were  somewhat  less  in 
terms  of  BCR,  the  trends  were  very  similar,  which  gives  one  confidence  in  the 
overall  validity  of  the  data. 

e.  Effect  of  Number  of  Layers  of  Reinforcement 

As  shown  in  Figure  53,  the  6CP  test  conducted  with  two  layers  of  woven 
fabric  reinforcement  indicated  somewhat  greater  BCRs  with  the  spacing  and  depth 
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less  than  1  inch,  whereas  the  increase  in  BCR  was  about  the  same  for  the  fabric 
layers  at  4  centimeters  or  1.6  inches  depth  and  spacing.  We  conclude,  there¬ 
fore,  that  the  benefit  from  two  fabric  layers  is  much  greater  if  the  depth  and 
spacing  of  the  reinforcement  layers  are  relatively  small  compared  to  the  diame¬ 
ter  of  the  loaded  area.  The  work  by  Binquet  and  Lee  (References  21  and  22)  sug¬ 
gests  that  the  failure  mode  that  would  be  observed  because  of  the  dimensions  of 
the  loaded  area  and  the  test  box  would  be  "ties  pulling  out."  With  short  ties, 
i.e.,  a  limited  distance  to  the  edge  of  the  reinforcement,  the  ties  do  not 
mobilize  sufficient  frictional  resistance  to  actually  exceed  their  breaking 
strength  in  tension;  rather  the  ties  will  simply  pull  out  of  the  sand  bed.  This 
is  the  type  of  behavior  that  we  observed,  especially  when  we  continued  to  push 
the  loaded  plate  into  the  sand  after  the  initial  bearing  capacity  failure  was 
observed.  For  evidence  of  the  reinforcement  pulling  out,  see  Figure  41. 

f.  Function  of  Type  of  Material 

The  tests  which  contrasted  the  woven  fabric  with  the  geogrid  reinforce¬ 
ment  were  the  3CP  and  PS  tests.  For  the  circular  plate  test,  the  BCR  for  the 
grids  was  about  3  percent  greater  than  the  fabric  BCR  at  d  =  1  inch,  and  about 
13  percent  greater  for  <1=2  inches.  For  the  plane  strain  test,  the  reverse  was 
true;  that  is,  the  grids  actually  had  a  lower  bearing  capacity  than  their  woven 
counterparts  at  both  depths,  although  the  differences  shown  in  Table  6  are 
small.  For  all  practical  purposes,  the  BCRs  in  plane  strain  for  both  types  of 
reinforcement  are  the  same. 

g.  tffect  of  Partial  and  Full-Edge  Fixity  of  the 

Rei nforcement 

As  shown  in  Figure  53,  edge  fixity  actually  decreased  the  BCR  for  the 
two  tests  conducted  at  d  =  1  inch.  This  response  is  surprising  because  it  is 
difficult  to  imagine  how  the  edge  conditions  would  influence  the  initial  bearing 
capacity  failure. 

Present  results  will  be  contrasted  with  those  obtained  by  Haliburton,  et 
al.  (References  12  and  13)  later  in  this  section. 

h.  Scaling  Effects 

As  described  earlier,  one  of  the  objectives  of  this  research  was  to 
determine  what  the  bearing  capacity  of  fabric-rei nforced  sands  would  be  under 
typical  fighter  aircraft  loads  and  tire  print  sizes.  For  reasons  explained  pre¬ 
viously,  it  was  not  practical  to  test  circular  plates  of  the  same  size  as  an 
aircraft  tire  print.  It  is  interesting,  therefore,  to  look  at  the  tests  in  an 
attempt  to  "scale  upwards"  to  the  approximate  9-inch  diameter  print  size  of 
typical  fighter  aircraft  tires. 

Figure  54  shows  the  maximum  load  on  the  plate  at  failure  vs.  the  ratio  of 
the  depth  of  burial  of  the  reinforcement  to  the  diameter  of  the  loaded  plate,  or 
d/D.  In  the  case  of  the  plane  strain  tests,  the  ratio  of  the  depth  of  burial  to 
the  width  of  the  plate  was  used.  For  the  sake  of  comparison,  the  theoretical 
Dearing  capacity  values  are  indicated  along  the  ordinate  by  small  arrows.  It 
will  be  recalled  that  these  are  the  theoretical  values  calculated  for  $  =  42 
degrees.  As  in  Figure  52,  there  is  a  marked  increase  in  maximum  load  due  to  the 
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presence  of  the  reinforcement.  In  every  case  of  d/D,  the  presence  of  the  rein¬ 
forcement  results  in  a  greater  ultimate  load  than  for  the  unreinforced  sand 
alone  for  the  same  size  plate.  For  the  case  of  the  plane  strain  tests,  the 
degree  of  improvement  is  significant.  In  Figure  55,  the  stress  at  maximum  load, 
°max/Aplate  ~  d»  is  shown  vs.  d/D.  The  trends  described  with  respect  to  Figure 
54  are  apparent  here  also. 

For  comparison  purposes,  data  were  extracted  from  Mali  burton,  et  al . 
(References  12  and  13)  and  their  data  (for  a  6-inch  plate  only)  are  shown  in  Fig¬ 
ure  55.  Shown  on  the  ordinate  is  a  small  arrow  labeled  "Th.B.C."  which  indi¬ 
cates,  as  before,  the  theoretical  (unreinforced)  bearing  capacity.  This  value 
is  for  y  =  35  degrees,  and  it  agrees  very  well  with  results  in  Reference  13. 
Actually,  the  trend  towards  increase  in  ultimate  stress  on  the  plate  due  to 
rei nforcement  is  about  the  same  as  found  in  the  present  study,  although  there 
appears  to  be  more  scatter  in  their  data. 

In  Figure  56,  the  BCR  is  shown  vs.  the  normalized  depth  of  the  reinforce¬ 
ment,  d/D.  Obviously,  for  no  reinforcement,  the  BCR  is  1.0,  so  the  points  for 
d/D  =  u  are  not  shown.  The  trends  as  suggested  by  Figures  54  and  55  are  also 
indicated  here.  It  was  expected  that  the  size  effect  for  the  circular  plate 
tests  would  be  small,  but  this  is  not  the  case,  as  shown  in  Figure  56.  The  3CP 
tests  have  a  higher  ultimate  BCR  than  the  6CP  test,  and  from  both  tests,  the  BCR 
decreases  as  d/D  increases.  As  noted  previously,  the  plane  strain  test  behavior 
is  in  the  opposite  direction.  The  responses  for  yrids  as  well  as  multiple 
layers  are  also  indicated  here.  With  Figure  56,  the  benefit  of  the  various 
reinforcement  systems  can  be  easily  ascertained  as  a  function  of  the  ratio  d/D. 

A  significant  effect  of  the  size  of  the  plate  on  the  BCR  has  been  observed 
by  others.  For  example,  Jarrett  (Reference  55)  conducted  several  large-scale 
laboratory  plate  bearing  tests  on  peats,  and  he  found  a  significant  difference 
between  tests  using  a  6-inch  diameter  plate  and  those  using  a  12-inch  diameter 
olate.  Thus,  using  purely  geometric  scaling  to  predict  the  load-carrying  capa¬ 
city  of  geotextile-reinforced  soils  is  not  possible.  Apparently  other  factors 
affect  the  results. 

For  comparison,  the  results  from  Haliburton,  et  al .  (References  12  and 
13)  are  shown  in  Figure  57.  Several  interesting  phenomena  are  apparent.  For 
one  thing,  there  is  a  rather  large  scatter  in  the  BCR  for  tests  at  d/D  =  0.5. 

For  the  tests  in  which  the  BCR  is  greater  than  1,  those  values  are  of  the  same 
order  of  magnitude  as  found  in  the  present  tests.  When  d/D  increases  to  1.0, 
there  is  actually  a  loss  in  efficiency  due  to  the  presence  of  the  reinforcement, 
and  this  is  true  also  for  the  4-inch  plates  at  d/D  =  0.5.  Why  this  effect 
occurs  is  not  obvious.  If  the  behavior  hypothesis  proposed  by  Haliburton  and 
his  colleagues  is  correct  (Section  I),  then  the  deeper  the  initial  level  of  the 
reinforcement,  the  less  effective  it  should  be.  In  other  words,  if  the  rein¬ 
forcement  is  initially  at  a  distance  b  or  greater  below  the  bottom  of  the  foot¬ 
ing,  the  rei nforcement  has  essentially  no  effect  on  the  deformation  and  failure 
pattern  of  the  footing  until  the  footing  settles  sufficiently,  or  well  past  the 
"ultimate."  This  effect  is  seen  in  the  work-hardening  of  Haliburton's  load- 
settlement  curves  and  by  our  data  in  Figure  51.  Haliburton’s  tests  were  con¬ 
ducted  with  edge  fixity,  whereas  almost  all  of  the  present  tests  were  not.  How¬ 
ever,  Figure  56  shows  that  edge  fixity  was  not  a  significant  factor  in  increas¬ 
ing  the  BCR.  For  some  unknown  reason,  it  may  actually  cause  a  decrease  in  the 


BCR.  It  was  originally  anticipated  that  Haliburton's  tests  would  provide  an 
upper  bound  to  the  level  of  improvement  due  to  the  reinforcement.  From  the 
results  in  Figures  56  and  57,  an  upper  bound  is  not  apparent. 

i.  Box  Size  and  Edge  Effects 

In  the  early  stages  of  this  research,  the  question  «js  raised  aoout  the 
relationship  between  the  size  of  our  loaded  plate,  at  that  time  6  inches  in 
diameter,  to  the  width  of  the  test  box  itself.  For  those  tests,  the  ratio  was 
about  5.  Thus,  the  test  series  with  the  3-inch  diameter  plate  was  planned  to 
assist  in  scaling  to  larger  plate  sizes  and  to  investigate  possible  edge 
effects. 

During  most  of  the  tests,  a  mechanical  dial  indicator,  accurate  to  0.01 
millimeters,  was  placed  outside  the  box  with  its  stem  at  the  elevation  of  the 
rei nforcement.  During  loading,  the  deflection  of  the  wall  of  the  box  was  noted. 
In  the  case  of  the  6CP  and  PS  tests,  the  maximum  observed  deflection  was  between 
0.2  and  0.35  millimeters.  No  wall  movement  was  detected  in  any  of  the  3CP 
tests.  It  was  apparent  that  the  observed  movements  correlated  directly  with  the 
loading.  That  is,  as  the  load  was  increased,  the  deflection  would  increase 
until  the  maximum  value  was  reached  at  approximately  the  time  of  initial  bearing 
capacity  failure. 

Because  the  observed  movements  were  so  small,  it  is  felt  that  any  edge 
effects  could  be  considered  negligible. 

j.  Haliburton's  "Optimum  Depth"  Concept 

As  described  in  Section  I,  Haliburton  and  his  associates  (References  12 
and  13)  proposed  a  concept  of  "optimum  depth"  for  explaining  the  increase  in 
bearing  capacity  due  to  the  presence  of  the  rei nforcement.  They  suggested  that 
the  optimum  depth  of  the  fabric  should  be  at  D.5  B  tan  <p.  For  the  present 
study,  the  most  probable  <j  =  42  degrees,  so  the  "optimum  depth"  for  the  3CP 
tests  would  be  1.35  inches  and  for  the  6CP  tests  2.7  inches.  These  values 
correspond  to  a  d/D  ratio  of  0.45  for  each  case.  The  results  in  Figures  52 
through  56  indicate  that  the  "optimum  depth/’  if  it  exists,  is  less  than  this 
value,  perhaps  even  less  than  1/3  D.  In  a  practical  sense,  however,  there  are 
severe  limitations  to  the  minimum  depth  of  placement  of  the  reinforcement,  even 
if  it  could  be  shown  that  the  "optimum  depth"  is  1/3  D  or  less.  A  practical 
minimum  placement  depth  is  probably  about  1  inch  or  perhaps  deeper,  especially 
if  ordinary  construction  equipment  is  employed  for  the  placement  of  the  top  sand 
layer.  Considering  this  limitation,  then  the  "optimum  depth"  of  placement  that 
naliburton  found  is  not  unreasonable  for  typical  tire  print  sizes. 

k.  Deflection  Basin  Measurement  and  Rut 

Development 

It  was  of  some  interest  to  determine  the  shape  of  the  deflected  surface 
of  the  sand.  These  results  are  shown  in  Figures  58,  59  and  60,  for  the  3C.P,  6CP 
and  PS  tests,  respectively.  In  these  figures,  the  surface  deflection  in 
thousandths  of  an  inch  is  shown  vs.  the  distance  from  the  center  line  of  the 
plate.  On  the  top  of  each  figure,  small  arrows  A  through  E  are  indicated. 

These  positions  correspond  to  the  locations  of  the  UCDTs  described  earlier  (see 
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Figure  58.  Shape  of  Deflection  Basin  for  the  3CP  Tests. 
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Figure  60.  Shape 


of  Deflection  Basin  for  the  PS  Tests. 


Figure  26).  In  dll  cases,  the  deflections  were  determined  at  one-half  of  the 
load  on  the  plate  at  first  yieldiny.  Considering  the  rather  small  deflections 
measured  (in  many  cases,  only  a  few  thousands  of  an  inch),  the  data  are  remark¬ 
ably  consistent.  It  is  difficult,  however,  to  see  any  significant  effect  of  the 
different  reinforcement  configurations  other  than  that  there  is  some  improve¬ 
ment,  i.e.,  some  decrease  in  deformation,  because  of  the  presence  of  th  •■ein- 
forcement.  Apparently,  the  stress  field  is  changed  by  the  presence  of  the  rein¬ 
forcement  well  before  yielding  or  "failure"  occurs.  This  is  consistent  with  the 
strain  diffusion  hypothesis  presented  in  Section  11. 

Note  that  the  shapes  of  the  deflected  surfaces  are  similar  to  those  shown 
in  Figures  d  through  14  for  the  strain  diffusion  hypothesis  (Section  II).  How¬ 
ever,  deflection  profiles  in  those  figures  are  in  terms  of  relative  deflection 
(relative  to  the  deflection  of  the  rigid  plate).  Therefore,  an  attempt  was  made 
to  normalize  the  deflection  measurements  in  Figures  58  through  60.  These  data 
are  shown  in  Figures  61  and  6?.  They  correspond  to  the  3CP  and  PS  deflection 
data  given  previously.  (Because  of  the  difficulty  of  interpreting  the  normal¬ 
ized  3CP  deflection  results  for  Figure  61,  a  similar  normal ization  procedure  was 
not  attempted  for  the  bCP  tests.)  In  these  figures,  the  deflection  at  P  7 
was  normalized  for  the  various  rei nforcement  configurations  with  respect  to  the 
plate  settlement  at  maximum  load,  which  occurs  just  before  the  plate  plunges 
into  the  subgrade  sand.  The  data  are  remarkably  consistent  in  the  plane  strain 
tests  (Figure  62).  However,  this  was  not  the  case  for  the  3CP  tests. 

It  can  be  determined  from  Figure  14  that  the  maximum  theoretical  heave  is 
u.2  at  x/a  =  3.  Figure  61  shows  that  for  the  3-inch  diameter  plate,  a  relative 
deflection  (heave)  of  0.2  occurs  at  r/b  =  1.3,  which  is  equivalent  to  x/a  =  2.6, 
or  very  close  to  the  theoretical  maximum  at  x/a  =  3.  For  the  6-inch  diameter 
plate,  the  maximum  relative  deflection  (heave)  was  0.42,  and  it  occurred  at  x/a 
•  1.5.  The  difference  between  the  3CP  and  6CP  behavior  is  probably  related  to 
some  unknown  scaling  effect  of  plate  diameter,  as  described  earlier  in  this  sec¬ 
tion. 

As  the  operational  criteria  of  some  aircraft  specifies  ruts  to  be  less  than 
about  1  inch  in  depth,  it  was  of  interest  to  determine  what  the  load  would  be  at 

I  inch  depth  for  the  various  loading  plate  sizes  and  reinforcement  configura¬ 

tions.  In  some  cases,  the  nature  of  the  load-settlement  curves  and  the  limita¬ 
tions  of  the  instrumentation  made  it  impossible  to  determine  exactly  what  the 
load  was  at  a  settlement  of  1  inch.  But  a  few  general  conclusions  can  be  drawn 

from  the  data  that  do  exist  (see  Table  6). 

In  a  large  number  of  cases,  there  really  was  very  little  difference  between 
Prux  (load  at  initial  yield)  and  P  at  s  ■  1  inch.  This  observation  was  true 
independently  of  whether  the  fabric  or  grid  reinforcement  was  used.  However, 
this  was  not  true  for  the  cases  with  partial  and  full  fixity.  The  estimated  P 

at  s  -  1  inch  was  at  least  twice  that  of  F\tUX  .  These  observations  tend  to  ver¬ 

ify  the  observations  of  Haliburton,  et  a  1  ..'References  12  and  13).  They  found  a 
significant  "work  hardening"  of  the  load-deflection  relationship  after  substan¬ 
tial  (0.3  to  0.5  0)  sinkage  of  the  loading  plate  took  place.  If  the  fabric  is 
"fixed"  at  the  edge  and  placed  relatively  close  to  the  surface  of  the  sand,  then 
it  could  be  expected  that  significant  increase  in  bearing  capacity  and  decrease 
in  deformations  before  failure  would  occur  (see  Sections  1.3,  III  .  5 .  h ,  and 

I I  i  .5.  j  of  this  report) . 
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In  several  "reinforced"  tests,  loading  was  continued  well  past  the  initial 
failure  so  that  the  loading  plate  penetrated  2  inches  or  more  into  the  subgrade. 
Then  the  top  layer  of  sand  was  carefully  excavated  and  the  reinforcement 
removed.  Significant  depressions  in  the  exposed  sand  surface  were  observed 
(see,  for  example.  Figure  35).  In  contrast,  Haliburton,  et  al .  (References  13, 
p.  72)  reported  finding  no  such  permanent  deformation  below  the  level  of  the 
reinforcement.  The  difference  between  present  test  observations  and  those  by 
Haliburton  was  probably  the  result  of  the  edge  fixity  in  his  tests,  but  this  is 
only  speculative. 


SECTION  IV 


IMPLICATIONS  FOR  ALTERNATE  LAUNCH  AND  RECOVERY  SURFACES 

The  experimental  and  analytical  investigations  described  in  the  previous 
two  sections  have  rather  important  implications  for  the  design  and  construction 
of  alternate  launch  and  recovery  surfaces.  This  statement  is  made  in  spite  of 
the  fact  that  none  of  the  rei nforcement  configurations  actually  tested  wcul  :  be 
able  to  support  present  military  aircraft  on  unpaved  but  reinforced  surfaces 
constructed  of  dense  sands.  The  results,  however,  do  point  the  way  towards 
several  practical  techniques  which  could  significantly  enhance  the  bearing  capa¬ 
city  of  unpaved  surfaces  and  may  ultimately  provide  economic  ALRS  systems. 

With  the  exceptions  discussed  in  Section  III,  the  test  results  presented 
herein  tend  to  verify  the  conclusions  of  Haliburton,  et  al .  (Reference  13)  as  to 
the  effect  of  lateral  restraint  reinforcement  where  membrane-type  support  is  not 
expected.  However,  if  a  very  soft  subyrade  exists,  then  one  would  expect  the 
greatest  effect  of  the  membrane  to  be  realized.  This  was  shown  theoretically  by 
Bourdeau,  Harr  and  Holtz  (Reference  44;  see  Section  I  and  Appendix  A),  but  it 
could  not  be  verified  experimentally  in  the  present  series  of  tests  because  only 
dense  sands,  not  soft  clays,  were  tested. 

Tne  important  effect  of  geometry,  that  is,  the  depth  of  placement  o'  the 
f  i  r-s  t  layer  of  rei  nforcement,  as  well  as  the  effect  of  multiple  layers  of  rein¬ 
forcement,  was  shown.  However,  the  concept  of  a  "critical  depth"  of  reinforce¬ 
ment  at  0.5  B  tan  was  not  verified,  although  the  data  are  not  sufficiently 
precise  to  recommend  an  alternate  numerical  value  of  "critical  depth."  If  it 
exists,  it  probably  is  less  than  0.5  B  tan  $  or  about  1/3  3. 

The  requirement  of  sufficient  depth  of  cover  for  proper  fabric  anchorage  is 
apparently  not  as  critical  as  had  been  originally  thought.  From  the  tests 
reported  in  Section  III,  even  shallow  depths  of  cover  are  sufficient  to  develop 
fabric  frictional  resistance.  This  result  is  shown  even  with  the  "lower  bound" 
tests  with  no  edge  fixity  in  comparison  with  the  results  of  partial  and  full 
fixity  at  the  same  depth  of  fabric.  This  observation  has  important  practical 
implications.  In  the  field,  the  depth  of  cover  to  the  first  reinforcement  layer 
can  be  relatively  shallow,  depending  on  the  depth  required  to  develop  sufficient 
bearing  capacity  versus  the  practical  consi derati on  of  how  to  actually  construct 
a  thin  layer  of  sand  over  a  geotextile.  This  means  that  the  fabric  can  be 
placed  at  a  relatively  shallow  depth,  contingent  only  upon  practical  construc¬ 
tion  constraints,  and  not  contingent  upon  the  requirement  of  a  sufficient  depth 
of  burial  to  develop  fabric  frictional  resistance.  The  value  of  the  geotextile 
rei nforcement  is  apparently  much  greater  for  ALRS  installations,  which  are  only 
required  to  support  a  relatively  few  operations  before  excessive  surface  defor¬ 
mation  or  rutting  may  occur.  Excessive  surface  deformation  indicates  that  a 
bearing  capacity  "failure"  has  occurred  and  that  the  soil  exists  in  a  state  of 
plastic  equilibrium.  Thereafter,  there  will  be  little  or  no  increase  in  bearing 
capacity,  even  with  very  large  deformations. 

Based  on  the  results  presented  herein,  the  system  most  likely  to  be  suc¬ 
cessful  as  a  reinforcement  system  for  ALRS  would  be  multiple  layers  of  a  heavy 
geotextile  or  geogrid  with  the  depth  to  the  top  layer  as  small  as  can  be  practi¬ 
cally  constructed.  The  minimum  depth  should  be  on  the  order  of  I  inch  or  less 
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if  it  is  possible  to  construct  such  a  thin  layer,  although  scaling  upward  to 
aircraft  tire  sizes  may  alter  this  recommendation.  Construction  problems  will 
exist  also  with  the  requirement  for  multiple  layers  of  reinforcement.  Specific 
design  rules  for  ALRS  systems  as  well  as  specification  of  the  desired  fabric 
properties  are  not  possible  at  this  time  without  the  results  of  experiments  simi¬ 
lar  to  those  proposed  for  Phase  II  of  this  effort  (Section  V). 


SECTION  V 


OUTLINE  OF  PLANS  FOR  PHASE  II  RESEARCH 

It  was  originally  envisioned  that  the  laboratory  and  analytical  portions  of 
the  research  would  allow  a  program  of  full-scale  field  tests  to  be  planned  and 
carried  out  as  a  second  phase  of  the  research.  The  program  would  assess  the 
validity  of  the  analytical  relationships  developed  in  Phase  I  and  would  gather 
additional  experimental  evidence  on  the  performance  and  economics  of  selected 
soil  reinforcement  systems. 

1.  FIELD  TESTS 

Full-scale  model  tests  of  a  limited  number  of  fabric-rei nforced  runway  sys¬ 
tems  would  be  constructed.  These  test  sections  would  be  constructed  to  simulate 
actual  field  installation  procedures.  The  size  of  the  test  sections  would 
depend  on  the  place  where  the  tests  were  carried  out.  If  conducted  at  the  USAE 
Waterways  Experiment  Station,  Vicksburg,  Miss.,  then  test  sections  on  the  order 
of  20  feet  wide  oy  150  feet  long  are  possible,  as  indicated  by  Webster  (Refer¬ 
ences  36  and  37).  Thus,  five  different  test  sections,  each  30  feet  long,  could 
De  investigated.  At  WES  it  is  possible  to  have  test  sections  on  both  very  soft 
clay  (CBR  <  1)  subgrades  or  subgrades  of  loose  sands  and  silts.  It  is  recom¬ 
mended  that  a  sand  subgrade  be  utilized  initially.  One  section  would  be  unrein¬ 
forced  as  a  control,  and  the  other  sections  would  consist  of  different  types  and 
geometries  of  reinforcement.  The  test  sections  would  be  instrumented  for  vert¬ 
ical  and  norizontal  deformations.  I nductance-type  soil  strain  gages  would  be 
attached  to  the  rei nforcement  so  that  horizontal  deformations  could  be  deter¬ 
mined  (Reference  56).  Settlement  profiles  could  be  determined  either  by  cross- 
sectional  surveys  or  transverse  settlement  yages  (Reference  57).  The  response 
of  the  various  reinforced  systems  would  be  obtained  by  use  of  the  F-4  load  cart. 

The  cart  would  be  repeatedly  towed  over  the  control  and  reinforced  sections 
until  either  subgrade  failure  or  deep  rutting  of  the  surface  occurred.  Prior  to 
failure,  surface  deflections  would  be  obtained  by  the  traveling  laser  beam  sys¬ 
tem  developed  by  Elton  (Reference  58)  presently  at  WES.  In  this  manner,  the 
number  of  passes  of  the  F-4  load  cart  could  be  correlated  with  performance  of 
the  reinforced  sections.  The  observations  would  verify  predictions  of  prototype 
performance  under  full-scale  traffic  conditions  and  allow  empirical  corrections 
to  oe  made  to  design  procedures  developed  from  the  analytical  and  experimental 
work  conducted  in  Phase  I  of  the  research. 

2.  ADDITIONAL  LABORATORY  AND  ANALYTICAL  INSTIGATIONS 

A  limited  series  of  laboratory  investigations  should  be  carried  out  as  part 
of  the  Phase  II  research  program.  Additional  plate  load  experiments  should  be  car¬ 
ried  out  on  (a)  clay  subgrades  anci  (b)  on  the  geogrids  with  coarse  granular 
material  such  as  pea  gravel  utilized  for  the  top  surface  layer.  As  mentioned 
previously,  one  of  the  advantages  of  geogrid  reinforcement  is  in  their  ability 
to  interlock  coarser  granular  particles  to  provide  tensile  resistance,  in  addi¬ 
tion  to  increased  frictional  resistance.  Additional  mul tiple-layer  reinforce¬ 
ment  systems  should  also  be  investigated  to  supplement  previous  research.  Tests 
on  soft  clay  subgrades  would  enable  behavior  and  performance  to  be  extended  to  a 
wide  variety  of  natural  soil  types. 


3.  ADDITIONAL  ANALYTICAL  WORK 


The  potential  of  the  strain  diffusion  theoretical  development  described  in 
Section  II  has  barely  been  realized.  Work  on  this  aspect  was  begun  rather  late 
in  the  project  and, consequently, ic  could  not  be  completed  in  the  time  remaining. 
This  technique  is  not  only  applicable  to  ALRS  and  other  pavement  systems  but 
also  to  an  entire  class  of  foundation  (settlement)  problems.  The  possibility  of 
being  able  to  avoid  having  to  determine  or  estimate  traditional  constitutive 
relations  for  a  wide  variety  of  soils  is  extremely  attractive.  It  is  therefore 
recommended  that  additional  analytical  work  also  be  carried  out  as  part  of  Phase 
II. 


SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 


1.  CONCLUSIONS 

a.  Significant  increases  in  bearing  capacity  due  to  both  the  depth  of  the 
reinforcement  and  the  number  of  reinforcement  layers  were  observed.  However, 
the  improvement  effect  decreased  as  the  depth  of  the  first  layer  increased. 

b.  As  observed  by  others,  a  strain-hardening  effect  due  to  the  presence  of 
the  rei nforcement  was  also  observed  in  these  tests. 

c.  At  shallow  depths  of  reinforcement,  an  increase  in  the  modulus  of 
subgraoe  reaction  was  observed. 

d.  The  concept  of  a  "critical  depth,"  if  it  exists,  was  found  to  be  some¬ 
what  less  than  0.5  B  tan  4,;  it  is  probably  about  1/3  B  (or  1/3  D). 

e.  For  the  geometries  tested,  behavior  of  reinforcement  systems  in  terms 
of  maximum  load  on  the  plate  and  bearing  capacity  ratio  when  loaded  in  plane 
strain  was  the  opposite  of  that  observed  for  the  circular  plate  tests. 

f.  Based  on  a  limited  series  of  tests,  edge  fixity  conditions  are  not  con¬ 
sidered  to  be  important.  The  bearing  capacity  tended  to  decrease  for  these 
tests. 


g.  The  benefit  of  mul tiple-rei nforcement  layers  is  greater  if  the  depth 
and  spacing  are  relatively  small  with  respect  to  the  diameter  of  the  loaded 
area. 


h.  The  behavior  and  degree  of  improvement  using  geogrids  were  about  the 
same  as  with  the  woven  polyester,  even  though  the  grids  have  a  significantly 
greater  tensile  strength  and  modulus  than  the  geotextile.  This  difference  was 
probably  due  to  the  lack  of  any  contribution  of  "interlock"  with  the  sand  used 
in  the  tests. 

i.  Geometric  scaling  of  load,  bearing  capacity,  etc.  was  not  possible; 
apparently  factors  other  than  diameter  of  loaded  plate  influence  the  response  of 
geotextile-reinforced  systems. 

j.  Surface  depression  profiles  due  to  the  loaded  plate  could  be  predicted 
approximately  from  the  strain-diffusion  hypothesis  developed  in  Section  2. 

2.  RECOMMENDATIONS 

Recommendations  for  additional  research  in  earth  reinforcing  were  outlined 
in  Section  V.  Full-scale  reinforcing  systems  should  be  tested  under  prototype 
loading  conditions.  Construction  problems,  as  well  as  field  performance,  can  be 
modeled  adequately  to  effectively  permit  prediction  of  prototype  behavior. 


systems  can  be  finalized. 
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APPENDIX  B 


LISTING  OF  COMPUTER  PROGRAMS 


PROGRAM  MAimiNPUT,0UTPUT,TAPE5=INPUT'. TAPES  OUTPUT) 

c 

c  ***••••••••••••••«••••«••••••••••••••••••••••••••••••••••• 

c 

C  •  THIS  PROGRAM  EMPLOYS  A  2D  SYSTEM  OP  TWO  LAYERS  0 F  * 

C  *  SOIL  REINFORCED  BY  A  MEMBRANE  AT  THE  INTERFACE.  « 

C  *  WINKLER  MODEL  FOR  THE  LOWER  LAYER.  • 

C  *  TANGENTIAL  INTERACTION  IS  RIGID-PLASTIC  * 

C  *  ITERATIUE  SOLUTION  TO  THE  FINITE  DIFFERENCE  EQUATION  * 

C  *  DIMENSION  OF  UECTORS  *  300  • 

C  *  MODE  BATCH  • 

C  *  * 

C  *•**«•*»*•*****•****«****»*•«********«##»*•#•******«*•*•*•• 
C 

c 

COMMON  /BL 1  /EPS  ILO 1 .  EPS  IL02 .  TO.  DELTAT ,  I TMAX.  I  TER.  UM. 


*  EM.L.WR, I OUT ( 30 ) «  JOUT . XLA. A. TMAX. XTMAX 
COMMON  /SIMPS/FC300). NS. HS. FINT 

COMMON  /MAT1/XAI 300. 3) . XB(300)«U(300).NP 
COMMON  /MAT2/T ( 300 ) . TAU( 300 ) . LAMBDA. H . 

COMMON  /TG1/TAU1MC 300 ) . TAU2MC300 ) . TAU1 ( 300) . TAU2C 300 ) . X(300) 
COMMON  /TG2/S1ZC300), 522(300) 

C 

REAL  KS.OUI.L.LDIF.NU1. LAMBDA 
C 

C  INPUT  DATA 

C  . * 

C 

DO  100  1=1,30 
100  10UTI 1 1=0 

C 

URITEIE.901) 

301  FORMAT!/,/, IX. 130 (IN* ),/./, 45X. ISHPROGRAM  FABRIC4. 

*  /./. IX, 130(1H*)./) 

WRITEIB.305) 

90S  FORMAT ( /, /, IX. 37HC0NUENTIQN:  TENSILE  STRESSES  NEGATIUE) 

READCS. 1103)A 
READ(5. II03JHI 
1103  FQRMAKF10.3) 

READ(5. I103)NU1 
READ<5.  U03)KS 
LAMBDA=KS 
READ(5. II03)EN 
READC5. 1103)P 
READIS, 1 103)L 
READ(5. I103)FI 
READ<5. 1 103)F2 
READCS. 1103IGAMA 
READ<5, 1115)NSTEP 
U15  FORMAT  ( IB) 

C 

200  H=L/NSTEP 
C 

C  NUMBER  OF  NODES 

C 

NP=NSTEP+1 
READ(5, 1103)TO 
READ (5, 1103)DELTATO 
READ(5. 1214)EPSIL01 
READ<5, 1214 )EPSIL02 


READ'S. 1218)ITERHflX 
READ'S.  1218>  UMAX 
1214  FQRMAT'FIO.B) 

1218  FORMAT ( 13) 

C 

C  INTERMEDIATE  OUTPUT  COMMAND 
C 

READ<5, 1219)N0UT 

1219  FORMAT  C 12) 

IF  (NOUT.  N£.  0>CO  TO  210 
GO  TO  230 

210  DO  220  1=1. NOUT 
220  READ'S. 1219) IQUT(I) 

230  CONTINUE 

C 

C  OUTPUT 

c 

c 

WRITE' G. 2001 1 

2001  FORMAT'/.  IX. 21HDATA  FDR  COMPUTATION: , 2iH«»**«***********'**»'<0 
C 

URITE(5. c  005) 

2005  FORMAT(4SX. 12HCENERAL  DATA./) 

WRITE'S. 1 1 1 1 JHl.NUl. KS.EM.P.L. F1.F2, A 
1111  FORMAT'/. IX, 3HH1=. E12. 4.  IX. 1HM.9X 

1  • 4HNU1=. E12. 4. 3X 

2  . 3HKS=. E12. 4. 1X.BHKN.M-3. 3X 

3  . 3HEM=. E 12. 4, IX. 6HKN.M-1 , 3X 

4  .  2HP=.  E12. 4,  IX.  lrtl.  10X 

5  . IX. 2HL=. EJ2. 4. IX. 1HM. 10X 

6  . 3HF1  =  «  E 12. 4. 10X 

7  , 3HF2=. El 2. 4. 10X 

8  .2HA=,E12.4) 

WRITE'S. 1212IGAMA.NSTEP 

1212  FORMAT (IX. SHGAMA=. E12.4, IX. BHKN.M-3. 2X. EHN5TEP=. 14) 

URITE'G.2101) 

2101  FORMAT'/. 40X.34HDATA  OF  THE  ITERATION  COMPUTATIONS) 

UR ITE ( G. 2102 ) TO. DELTATO. EPSIL01 . EPSIL02.  ITERMAX.  ITMAX 

2102  FORMAT'/. IX. 3HTQ=, £12. 4, IX. GHKN.M-I . 4X, 8HDELTAT0=. El 2. 4. IX 

1  , GHKN.M-1,/, IX. 8HEPSIL01=. E12.4, BX. 8HEPSIL02=»E12.4. 5X 

2  , 8HI TERMAX®, 14. 13X. 6HITMAX*. 14) 

C 

C 

c 

C  INITIALIZATION  OF  UECTORS 

C  ••*•»»«•*••••*•••••»*♦•**• 

c 

DO  300  1=1.300 
TAU( I )=0. 

TAU1( I )=0 . 0 
TAU2C I )=0. 0 
TAUIM' I )=0.0 
TAU2M' I )=0. 0 
T( I )=0.0 
XA  ( 1 , 1 )  =0 . 0 
XA( I.2J-0.0 
XA(I.3)=0.0 
X8( I )=0. 0 
F( I )=0. 0 
X(I)=0.0 
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WCI)=0.0 
S22in=0.0 
300  S1Z( I ) =0. 0 

c 

C  CALCULATION  OF  TAUlfl 

c  . . •••••« 

c 

Q=P/(2«A> 

DO  3S0  1=1, NP 
X( I)=(I-1)*H 

Y1=(X(  I)-A)/(H1«S0RT(NU1)) 

Y2=( X( I  )+A)/(Hl»SQRT (NU1 ) ) 

HT=(Y2-Y11/NSTEP 
DO  320  J=1.NP 
F(J)=Y1+((J-1)«HT) 

F(J)=£XP(-(F(J)««2)/2) 

320  F(J)=F(J) /SORT ( 8* AT AN ( 1.0)) 

HS=HT 
NS=NSTEP 
CALL  SIMPSON 
350  S1Z( I )=FINT*0 

C 
C 
C 

CALCULATION  OF  DELTAP 
DO  370  1=1. NP 
370  F(I)=S1ZCI) 

NS=NSTEP 

HS=H 

C 

CALL  SIMPSON 
URITE(G.  12211FINT 

1221  FORMAT!/. IX. 7HPSICL)=.E12.4) 

DELTAP=P-(2«FINTJ 
URITECG. 12221DELTAP 

1222  FORMAT ( IX. 7HDELTAP=.E12.4, 7H  KN.H-1) 

DO  400  1=1. NP 

TAU1M(I)=<S1Z<I)-KGAMA«M1))*F1 
400  CONTINUE 
C 

C  UMAX  UITHOUT  MEMBRANE 
C 

m=sizcn  /LAMBDA 

c 

c  COUNTER  OF  PRINCIPAL  ITERATIONS  ITER 
C 

ITER=0 

JOUT=l 

C 

C  ITERATION  COMPUTATION  FOR  MAIN  CYCLE 
C 

450  ITER=ITER+1 

DELTAT=DELTATO 

C  _ 

C  FIRST  ITERATION  :  COMPUTATION  WITHOUT  TANGENTIAL  FORCES 
C 

IF  (ITER.  NE.  1)  GO  TO  530 
C 

500  DO  520  1=1. NP 

520  T(I)=TO 
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c 

C  1TERATIUE  COMPUTATION  FOR  SECONDARY  CYCLE 
C 

S30  CALL  ITERA 

C 

C 

C  CALCULATION  OF  S2Z  AND  TAU2N 
C 

DO  5S0  1=1. NP 
S2Z( I )=U(I) -LAMBDA 

5S0  TAU2M ( I  )=(S2Z(  I )  +  ( CAf1A«Hl )  ) »F2 

C 

C  COMPUTATION  OF  TANGENTIAL  FORCES 
C 

CALL  TANG 
C 

C  COMPUTATION  OF  EFFECT 1UE  LENGTH 
C 

IF  (FI.  NE.  0)  CO  TO  5G5 
IF  (F2.  NE.  0)  GO  TO  5G5 
XLA=L 
GO  TO  575 
565  CONTINUE 

C 

DO  570  1=1. NP 

IF  (TCI).  NE.  0)GO  TO  570 

XLA=X( I ) 

GO  TO  575 
570  CONTINUE 

C 

575  CONTINUE 

C 

C  CONTROL  OF  THE  MAIN  CGNUERGENCE  ON  UO 
C 

IF  (ITER.  EO.  1 )GO  TO  600 
DELTAU=(UOL-U(l))/U(l) 

ADELTAU=ABS( DELTAU) 

IF  (ADELTAU.  LT.  EPSIL011G0  TO  680 
C 

C  COMPUTATION  TOO  LONG 

c 

IF  (ITER.  EO.  ITERMAX)  GO  TO  900 
600  UOL=U( 1 ) 

C 

C 

C  INTERMEDIATE  OUTPUT 

c  . . 

c 

IF  (ITER.  NE.  TOUT ( JOUT ) )GO  TO  620 
CALL  OUTPT 
JOUT=JOUT«-l 
620  GO  TO  450 

C 

C  IF  CONUERGENCE  TAKES  PLACE  END  OF  COMPUTATION 
C 

680  CONTINUE 

C 

C  CALCULATION  OF  T(X) 

C 

N=NP-2 
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700 


DO  700  1=1. N 

F ( I )=( -3«W( I ) H(4«W( 1+1) )-H( 1+2) 
Fa)=S3RTU-KF(I)/(H»H))*«2) 

F(NP-l)  =  <U(NP)-il(NP-l)).'H 
F(NP-n=S0RT(l+F(NP-l)**2) 

F(NP)=r (NP-1) 

DO  720  1=1. NP 
720  T( I )=T ' I ) »F( I ) 

C 

C  CALCULATION  OF  TMAX 
C 

TMAX=T(1) 

XTMAX=0.0 
DO  730  I=2»NP 

IF  C T(  I ) .  CE.  TMAXIGO  TO  730 
TMAX=T ( 1 ) 

XTMAX=X(I) 

730  CONTINUE 

C 

UR=U(1)/UM 
COLL  OUTPT 
C 

GO  TO  393 
C 

C  COMPUTATION  TOO  LONG 
C 

900  WRITECG. 1801) 

1801  FORMATS. 1X.32HN0  CONUERGENCE.  IN  THE  MAIN  CYCLE) 

GO  TO  S80 

999  WRITE! G. 9001 ) 

9001  FORMAT (x. IX. 18HEND  OF  COMPUTATION) 

C 

STOP 
END 
C 

c 

SUBROUTINE  ITERfl 
C 
C 
C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C  COMMON  .'BH/EPSIL01,EPSIL02,T0»DELTAT.  UMAX,  ITER.UM.EM.L.WR 

C  1  I CUT ( 30 ) » JCUT • XLA. A. TMAX. XTMAX 
C  COMMON  /MAT1/XA(300>  3) . XBC300). U(300) »NP 

COMMON  /MAT2/T (300 ) , TAU( 300 ) , LAMBDA. H 
COMMON  ✓S IMPS'E ( 300 ) . NS. HS>  FINT 
COMMON  XTG2^S12(300).S22(300) 

COMMON  /TG1^TAU1M<300).TAU2M(300).TAU1(300).TAU2(300).X(300) 


SUBROUTINE  ITERA 

CYCLES  OF  SECONDARY  ITERATION  OF  FABRIC* 

SOLUTION  OF  FINITE  DIFFERENCE  EQUATIONS  On>E 
UERTICAL  AND  HORIZONTAL  EQUILIBRIUM 
FROM  THE  RESULTS  OF  THE  PREUIOUS  PRINCIPAL  ITERATION 
FORM  THE  MATRIX( SUBROUTINE  MATRIX) 

SOLUTION  OF  THE  LINEAR  SYSTEM(SUBRTN.  SLIMTRI) 
CONTROL  OF  THE  SECONDARY  CONUERGENCE 
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non  u)  nnono  non  a)  norm  onn  —  nnno  n  non  oonon  no 


T 

-  V-*-  ■ 


PEAL  LAMBDA. LD1FF.L 
INITIALIZE 

COUNTER  OF  THE  SECONDARY  ITERATIONS  IT 
I  T  =  1 

COUNTER  OF  THE  JUMPS  OF  THE  ROOTS  DURING  CONUERGENCE 
XSAUT  =  0. 0 
DIFF-Q . 0 

DIFFERENCE  FROM  THE  PREUIOUS  ITERATION 
00  LD1FF -DIFF 

FORMAT  ION  OF  THE  MATRIX  COEFFICIENTS 
CALL  MATRIX 


SECOND  MEMBER  OF  THE  SYSTEM 

NSTEP=NP-1 
DO  200  I =2. NSTEP 
oo  xb;I)=siz(I)»h«h 
XB( 1 ) =0 . 0 
X3(NP)=0.0 

SOLUTION  OF  THE  SYSTEM 
CALL  SLIMTRI 
CONTROL  OF  CONUFPGENCE 
CALCULATION  OF  SRA 
N=NP-2 

DO  300  1=1. N 

F(I)  =  (-3«L'(I))*(4*U(I  +  l))-U(I+2) 

oo  F(i)=soRT(i  +  (F(D/(a*Hi)**a) 

F (HP-1  )  =  (UffiP)-U(HP-l  >  )/H 
F(NP-1  1 =SORT  < 1 *F l  MP-1 )*»2) 
F(NP)=F(NP-1 ) 

NS=NSTEP 

H5=H 

CALL  SIMPSON 
SRA=F 1NT-L 

CALCULATION  OF  SRB 

DO  350  1=1. NP 

350  F<  I  )=F(  n»ABS(T(I) ) 

NS=N5TEP 

HS=H 
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CALL  SIMPSON 
SRB=FINT/EM 
C 

C  TEST  OF  CONUERGENCE 
C 

DIFF=SRA-SRB 

ADIFF=A5S(DIFF) 

ADIFF  =ADIFF/SRA 

IF  (ADIFF.  LT.  EPSIL02)G0  TO  BOO 
C 

C  TESTS  ON  THE  JUMPS  OF  THE  ROOTS 
C 

PROD=DIFF*l.DIFF 
IF  (PROD.  LT.  0)  GO  TO  550 
C 

C  STEPS  OF  THE  JUMP 
C 

540  CONTINUE 

TO=TO+DELTAT 
DO  54S  1=1. NP 
545  T(  I )=T( I )+DELTAT 

I T  =  1 T  + 1 

IF  (IT.  GT.  ITMAX1GO  TO  580 
GO  TO  100 
C 

C  JUMP 
C 

550  XSAUT=XSAUT*1 

XSIGNs-1 

570  D£LTAT--XSICN»(C1-'(2»XSAUT))«DELTAT) 

TO=TO+DELTAT 
DO  575  1=1. NP 
575  T( I )=T( I l+DELTAT 

IT=IT+1 

IF  (IT.  GT.  ITMAXJGO  TO  580 
C 

C  INTERMEDIATE  OUTPUT 

c 

IF  (ITER.  NE.  IOUT( JOUT) )GO  TO  578 
WRITEfS.BOODITER 
2001  rCRMAT(x, 1X.SHITER=.I4) 

•TE(G. 1001 IDELTAT 
1001  .  'MAT (IX. 7HDELTAT=.E12.4) 

URITEtE. 1411) 

1411  FQRMAT85X. 4HSRA=, E 12. 4, 5X, 4HSRB=.E12. 4.5X. 5HDIFF*. 

1  E12.4.5X. EH AD I FF=. El 2. 4) 

URITE(G. 10 12) PROD 
1012  FORMAT (IX. 5HPROD=. E24 .4) 

578  GO  TO  100 

C 

580  URITECG. 1420) ITER. IT 

1420  FORMAT (/.  /.  IX.  37HNO  CONUERGENCE  IN  Tr.E  SECONDARY  CYCLE.  IX. 

1  5HITER=, I4,3HIT=I4) 

STOP 

C 

600  CONTINUE 

RETURN 
END 
C 
C 


l 
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SUBROUTINE  TANG 


MimmimWIWWHtHHMWHHHWHWIHUlHIWIH 

SUBROUTINE  TANG  • 

COMPUTATION  OF  THE  NODAL  TANGENTIAL  FORCES  FOR  * 

FABRIC4  * 

PRIMARY  TANGENTIAL  STRESSES  SIZ  AND  S2Z  * 

SECONDARY  TANGENTIAL  STRESSES (ELASTO-PLASTIC)  * 
TENSILE  NODAL  FORCES  IN  THE  MEMBRANE  T  • 


COMMON  /  BL 1  /  EPS  1 1.01,  EPSILC2.  TO.  DEL  TAT .  ITMAX.  ITER.  HM.  EM.  L,  HR 
1  .  I  OUT  ( 30 ) .  JC-UT .  X'l_  A.  A.  TMAX.XTMAX 

COMMON  /TGI / TAU1MI 300) . TAU2MI 300) . TAU1 (3001 . TAU2I300 ) . XC300) 
COMMON  / TC2/S1 Z ( 300) . S22(300) 

COMMON  /f1AT[/RA(  300. 3).  XB<  300). Ui  300) .IIP 
COMMON  /MAT2/ T (3C0).TAU( 300). LAMBDA, 

COMMON  /SIMPS/F( 300). NS. HS. PINT 

REAL  LAMBDA 

SECONDARY  FORCES  ARE  ZERO  AT  THE  EXTREMES 


F ( 1 )=0.0 
TAU1 ( 1 )=0. 
TAU2(l)-0. 
TAU(  n=o. 

F ( NP ) =0 . 
TAUCNP'/=0. 
TAU1 ( NP ) =0. 
TAU2(NP)=0. 

FINT=0. 

TINT^O. 

NSTCP=NP~l 


DO  500  I-2.NSTEP- 

CCMPUTATION  STARTS  PROM  THE  END  OF  THE  MEMBRANE 

AUXILL IARY  INDEX  IA 

I A=NSTEP-I +2 
JA=IAU 

LIMIT  OF  THE  MOBILIZED  STRESS 

RIGID  PLASTIC  MODEL 

•ISO  TAUHIA)=TAU1M(IA) 

TAU2(IA)=TAU2M(IA) 

500  TAU(IA)=TAU1(IA)+TAU2(IA) 

C 


J 


c  TENSILE  FORCES  IN  THE  MEMBRANE 
C 

DO  SOD  I=2*NP 

B=(TAU<  Jl+TAUCI)  )*H/2 
TINT =T INT  +  B 
t  a  )=th)»tint 

r 

C  NO  COMPRESSIONS 

"  if  am.  a.  qigo  to  eoo 

T( I  )  =  0. 

Tf)U(  I  )=0. 

TPU2< I )=0. 

600  CONTINUE 


PETUPN 

END 


r 

C 

c 


:<<■><  1. 1  -'0. 
W(1.2!--l. 
MAC  I • 3'-~l • 
y;M(NP.  1  i=0. 
Xh;NP,2)=1  . 
..h(NP,3)  0. 


Rc.TURN 

END 
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noon  nnonnnonnnn  no  nru  non—  noon  nnnoonnnonoo 


c 

SUBROUTINE  SLINTRI 


*  4 

*  SUBROUTINE  SLIMTRI 

*  SOLUTION  BY  THE  DIHINISHON  OF  ft  TRI-DIAGONAL  LINEAR 

*  SYSTEM  ' 


COMMON  /MATUXAtSOO.  3) . XBC300 > .  UC  300) .  NP 
ELIMINATION 


DO  100  1=2. NP 
FACT=XA< I. 1 )/XA( 1-1.2) 

XA( I . 2 ) =XA( I,2)  +  (XA(I~1. 3)*FACT) 

00  XB( I )=XB( I )  —  C  XB  C I— 1 ) •FACT ) 

COMPUTATION  OF  UC I )  BY  SUBSTITUTION 

U(NP)=XB(NPVXA(NP,2) 

□0  200  J=2. NP 
I=NP+1-J 

u<  n  =  <xB(  n-(XA(i.3)*ua+i)))/'XACi,2) 

RETURN 

END 


SUBROUTINE  SIMPSON 


•  SUBROUTINE  SIMPSON  * 

*  NUMERICAL  INTEGRATION  USING  5IMPS0N  RULE  • 


COMMON  /'SIMPS/F  (300) .  NS.  HS.  FINT 
SIMPSONS  FORMULA  (NS  HAS  TO  BE  EUEN) 


100  FINT=0. 

DO  110  J=2. NS. 2 
XINT=F( J-1)+4«F(J)+F( J+I) 
XINT=XINT«HS^3 
110  FINT=FINT+XINT 

C 


RETURN 

END 


nnnnnnnnnn 


SUBROUTINE  OUTPT 

. . . . 

SUBROUTINE  OUTPT  * 

OUTPUT  FOR  FABRIC4  * 

»»»•••!  •»••••*•••*••«-•••«*•**«•**•**•*•************”**** 


c 


COMMON  /BLl/'EPSILOl.EPSILOS.TO.DELTAT,  ITMAX.  ITER.  UM,  EM.L.UR. 

.  IOUTC30).  JOUT.>3_A,A.TMAX.XTMAX  _ftn.  ynnn. 

COMMON  ✓  TG1./TAU1H(300).  TAU2M1 300) » TAU1 ( 300 ) » TAU2(300) » XC300 ) 
COMMON  /TG2/S12<300). 522(300) 

COMMON  /MAT1''XA(300.3).XB(300).U(300).NP 
COMMON  /MAT2'T ( 300 ) ■ TAU( 300 ) » LAMBDA. H • 


PFAL  LAMBDA. L 


C 

1515 


C 

C 

1525 


300 

1526 

C 

C 


UR  I TE  r  6 . 1515)1 TER. U(1).T(1).T( ^  *  ™AX.  XTNAX.  XLA,  J® 
FORMAT ( ✓.  /. IX. 27HRESULTS  OF  THE  COMPUTATION. . 5X. 5HITER  » IS.- 


1  ✓  . 

2  /.25X.SHUC 1)=.E1P.4.SX. 1HM. 

3  25X.  5HT (  1 )=. E12.4. SX.  6HKN.M-1* 

4  /.25X.5HT(Ll=.E12.4.5X. 6HKN. M“1 » 

5  /. 25X. 5HTMAX=. E12. 4. 5X. GHKN.tf-l. 

6  /.<?5X.6HXTMAX=,E12.4.5X.  1HM. 

7  / ,  25X. 4HXLA-. E 12.4. 5X. 1HM. 

8  /.25X. 3HUM-.E12.4. 

9  /,25X.3HUR=,F6.3.) 


F0RnATM»?lHN.7X.  1HX.  10X.3HS1Z.  1  OX, 3H52Z. 3X. 4HTAU1. 9X. 
1  4PTAU2. 10X.3HTAU. 12X. 1HT. 12X* 1HU.X) 

UPITEtS. 1526) I.Xd ).S12(I).S2Z( I)»TAU1( I)TAU2(I)f TAUCI) 

1  . r(I).U(I)  „ 

FORMAT ( IX. 13. 1X.E12.4.7(1X.E12.4)) 


RETUPN 

END 


no  on  nnn  no  on  n  nn  nru  nru  *—  nnnor 


PROGRAM  STRAINC  INPUT,  OUTPUT.  TAPE5=  INPUT,  TAPES=OUTPUT> 

C 

c 

c  . . . . . . . . 

c  •  • 

C  •  PROGRAM  STRAIN  • 

C  •  • 

C  •  THIS  PROGRAM  COMPUTES  THE  SURFACE  SETTLEMENT  UNDER  • 

C  •  AN  APPLIED  CIRCULAR  STRAIN  IMPRESSION.  USING  STRAIN  • 

C  •  DIFFUSION  TECHNIQUES.  • 

C  •  • 

C  . . . . . . . . . . 

DIMENSION  C ( 250  )• DI (40 1 
REAL  MU 
C 

PI=3. 1427 

Q=  1.00,  ! ■  ’it,  i  rv  'Wi  ■  :  u  ,1  ■  i  ,  |  Hj  ■  •  i  j  i  ;  i  |  ii  ,  ;  i  ,  [  i  , 

’ t  /  :■  'in  i  -  ns.  'i.  :  \ •  v  Ai  .  n  i «  i  ».  >.i"  I.  iMA :  r'i  ■ 
mi"!-  mi  i  ,  ii  \  !f'iM  i)i  i  !■'  ii '  ; 

INPUT  MU  UALUE 

READ  (S.IO.'HJ,  ■ 

0  FORMAT  CFIO.3) 

URI TP(G.  20 1  MU 

0  FORMAT (•''/'» IX, 8H  MU  UALUE.F5.3) 

URITECG.25) 

5  FORMAT (  ", 6X, 3HR/A, 13X. 1GHSETTLEMENT  RAT TOT 

DO  600  M=l.  1G 
R  =  1.0  +  aW)*0.25 


R0=0.0 

DO  GOO  J=l,  •, 


Z=(J-1)«DEL2 
E=0. 5»  MU*(Z)*«2 
D=0.0 

DO  300  1=1.240 
DELX=0. 125 
X=(I-U»DELX 
A2=EXP(-(X**2)*B) 


IF  (X.  GT.  7 . 0)GO  TO  270 

USE  BESSEL  FUNCTION  FOR  LOWER  ARGUMENTS 

Sl=1.0 
S2=l .  0 


DO  200  K=1 • 10 
Nl  =  l 


CALCULATE  THE  FACTORIAL 


7 


non 


c 


DO  100  L=1.K 

N1=L»N1 

CONTINUE 


100 
c 
c 

Nii=(Ni)«*a 

c 

c 

XI -R*X 

DEL5=((-0.25»(X*»2))**K VNll 
DELS1=(  <-0.25*(Xl*»2))**K)'Nll 
C 

S1=S1+DELS1 

C 

DE1.S2=D£LS/(K+1) 

C 

S2-S2+DELS2 

C  IE  (ABS(DELSl).  GT.  0.0001)  GO  TO  150 

IF  (ABSCDELS2).  GT.  0.0001)  GO  TO  150 
C 

GO  TO  250 
C 

ISO  CONTINUE 

C 

200  CONTINUE 
250  CONTINUE 

C 

S2=S2»0 . 5*X 
C 

IF  (XI.  LT.  7.0)GO  TO  290 
C 

270  CONTINUE 

USE  APPROX I NOTE  BESSEL  FUNCTION  FOR  LARGE  ARGUMENTS 

THETA=X-PI'4.0 

S1  =  (S0RTC2/PI/X)  )*COS( THETA") 

IF  (X.  GT.  7.0)  GO  TO  290 
C 

THETA1=X-3.0*PI/4.Q 

S2= ( SORT ( 2/P I /X ) ) *COS ( THETA J ) 

290  CONTINUE 

C 

C 

C(I)=(A2)*(SI)*(S2) 

IF  (I.  EQ.  1 )GO  TO  300 
C 

C  C3HPUTE  THE  STRAIN 
C 

DFLD=0.5»(CCI)+CCI-1))*DELX 

D=0ELD+D 

IF  (ABS(DELD).  LT.  0.0001)  GO  TO  400 
C 

300  CONTINUE 

C 

400  CONTINUE 


C 

D1(J)=Q*D 

c 

IF  (J.  EO.  l)GO  TO  BOO 
C 

C  CALCULATE  THE  SETTLEMENT 

c 

DELf?0=0.5*(Dl(J)+Dl(J-l))*DELZ 

RO=RO+DELRO 

IF  (ABS(DELRO).  LT.  0.0001)  GO  TO  700 
C 

600  CONTINUE 
C 

700  CONTINUE 

C 

C  CALCULATE  THE  SETTLEMENT  RATIO 
C 

IF  (R.  EO.  1 . 000 )R01=RO 
R0=R0/R01 
C 

URITECG.750)R.RO 
750  FORMAT (/«F10.3.9X>F10.3) 

C 

C 

800  CONTINUE 

C 

STOP 

END 


1  * 

(Tile  reverse  of  this  paqe  is  blank) 


DATE 

FILMED 


